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Abstract 



A fundamental prediction of the Dirac theory for electrons in strong electric fields 
is examined. When the electron's binding energy is sufficiently large the state 
becomes a pair creation resonance state, which is called a supercritical resonance. 

This resonance is similar in character to other resonances in atomic physics. It is 
parameterized by its energy and its lifetime. A numerical discretization technique, 
the mapped Fourier grid method (MFG) is extended to the Dirac equation and is 
used to solve for the resonance parameters of a quasimolecular supercritical IScr 
state which arises, e.g., in a uranium- uranium collision. Direct methods using 
only the MFG method are shown to give reasonable estimates for the resonance 
parameters. Analytic continuation methods such as complex scaling (CS) of the 
coordinate or adding a complex absorbing potential (CAP) are then applied. They 
allow for more accurate determinations of the supercritical resonance parameters. 

The analytic continuation methods of smooth exterior scaling (SES) and CAP 
are then augmented by a Pade approximant extrapolation. Both extrapolated 
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methods are shown to give highly accurate results, with substantially smaller ma- 
trix sizes than direct analytic continuation calculations. They are also shown to 
give results in close agreement with each other. 

The (extrapolated) augmented analytic continuation methods are used to inves- 
tigate the effects of higher-order couplings, beyond the monopole approximation. 
For the nearly charge symmetric system of U^^"''-Cf^^+, it is shown that the S-D 
quadrupole coupling is the next dominant interaction after the monopole interac- 
tion. 

Collisions near the Coulomb barrier, in which supercritical resonance states 
occur, are also calculated using the MFG. Results are presented for the collision 
of U^^"'"-U^^^ for a zero-impact-parameter Coulomb trajectory at a center-of-mass 
energy of 740MeV. The enhancement effects to the positron spectrum due to nuclear 
sticking at closest approach is examined. Results are given for sticking times of 
2,5,10x 10~^^s. Nuclear sticking is shown to offer the possibility of demonstrating 
experimentally the existence of supercritical resonance states. 
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1 Introduction 



In nature, two types of particles are found and are distinguished by their multi- 
particle statistical behavior. Bosons are particles obeying Bose-Einstein statistics, 
and allow any number of particles to be in the same quantum state. Fermions 
obey Fermi-Dirac statistics. Fermions must all be in a unique quantum state and 
this is responsible for the structure observed in atoms, nuclei and in baryons. The 
fundamental building blocks in nature are fermions, namely the quarks (the neu- 
tron, proton and other baryonic matter are composed of quarks and anti-quarks) 
and electrons (including muons and tauons). The fundamental bosons appear as 
mediators of interactions (photons and gluons), but composite bosons, made up of 
an even number of fermions, also exist (e.g. the deuterium nucleus, "^H). While 
atomic physics is concerned with modeling the behavior of all atomic systems, the 
electron plays a central role, since it is stable, abundant and the lightest charged 
fermion making it very accessible experimentally. 

The electron has an intrinsic spin of | i.e., a half- integer value as do all fermions. 
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Bosons possess integer spin (e.g. photons have an intrinsic spin of 1). A notable 
characteristic of particles is that different particles with the same spin behave in 
the same manner by obeying the same quantum mechanical equations of motion. 
For instance, the muon is also a spin-i particle with the same electric charge as the 
electron, but it is more massive (around 207 times the mass of the electron) and 
has a finite lifetime (primary decay channel: /i~ — > e~z7gZ/^). Prior to decaying, the 
muon will behave as an electron with a heavier mass. Therefore, any model of the 
electron will also describe other fundamental (point) particles with spin-|, once the 
constants such as the mass or sign of the charge are adjusted. This includes the 
anti-particles of the fermions. The positron is the anti-particle of the electron. It 
has the same mass and spin, but the opposite charge. 

In order to describe the behavior of the electron it is necessary to specify its 
equation of motion. Initially the electron was modeled as a classical charged point 
particle with quantized orbital angular momentum. This was hypothesized by Bohr 
to explain the level structure of the hydrogen atom which consists of an electron 
orbiting a proton. The Bohr model of the atom was relatively successful and ex- 
plained the Balmer series (spectral lines of hydrogen), but fell short in many areas, 
e.g. when applied to the spectrum of helium or other many-electron atoms. The 
Bohr-Sommerfeld model, incorporated elliptical orbits, and addressed the radiative 
decay problem, but could not resolve the helium problem. It was replaced by the 
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model of an electron as a wave obeying the Schrodinger equation. Pauli augmented 
this theory by including the intrinsic spin of the electron, ad hoc, in what is now 
called the Schrodinger-Pauli equation. This is a useful model of the electron that 
describes most the quantum phenomena dealing with the electron. 

When describing the electron as a wave, each orbital has room for two electrons 
of opposite spin due to the spin differentiating the quantum state. Since the electron 
is a fermion, once an orbital is occupied by a pair of electrons no other electron 
can be in that state. This forces other electrons into higher states, and causes 
higher orbitals to be filled giving structure to the atomic system. The blocking of 
an orbital, due to being filled with an electron, is called Pauli-blocking. 

For all the success of the Schrodinger-Pauli equation, it was not the final model 
of the electron. This equation accounts for the fact that the electron is a fermion 
but lacks the effects of special relativity. This equation is, therefore, only valid at 
low energies {E ^ mpC^ where mg is the mass of the electron and c is the speed 
of light in vacuum). A model which takes into account special relativity must be 
valid for kinetic energies that are comparable to, or larger than, the rest energy of 
the electron. With such sizable kinetic energy it is possible to create an electron- 
positron pair {E ^ 2mcC^). A relativistic theory of the electron, which would to 
fulfill these requirements, is needed. 



3 



1.1 Relativistic quantum mechanics 

To satisfy special relativity, the total energy of a free particle must be equal to 
sjp^c^ + m^c^. It is also necessary to have the same-order derivatives for the coor- 
dinate and time dependence. Dirac demanded that they be of first order to have 
simple time evolution. This resulted in a multi-component wave equation, with four 
components being the lowest number able to satisfy Poincare invariance. This was 
the start of relativistic quantum mechanics (RQM) for fermions. RQM describes 
the electron at very high and very low energies. Using hole theory, RQM also ac- 
counts for the possibility of electron-positron pair creation. To accomplish this, the 
tools of quantum field theory (QFT) are used. 

In QFT, the classical field degrees of freedom are quantized. This is done by ex- 
pressing the classical field in normal modes and describing each mode by a harmonic 
oscillator equation. Canonical quantization is then used to promote the expansion 
coefficients to creation and annihilation operators. The first complete quantum 
field theory was created, by Dirac, in 1927 and described the electromagnetic field. 
By combining the Dirac equation for the electron with the Maxwell equations for 
the photon and quantizing them, by imposing anticommutation and commutation 
relations respectively on the field operators, a robust part of the model of charged 
fermions is obtained. This is called quantum electrodynamics (QED). In QED the 
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photon acts as the force carrier for the electric charge, communicating the interac- 
tion between the different particles. 

RQM predicts that an electron-positron pair can be created both dynamically, 
through excitation by a time-varying field, or statically though very strong poten- 
tials. While the former is well understood and verified by experiment, the static 
creation of pairs has yet to be demonstrated experimentally. This method of pair 
creation is intimately tied to the fundamentals of QED. Thus, it presents a good 
test of the theory in the presence of very strong potentials. 

In RQM, the spin-^ particles are described by the Dirac equation. This is done 
in the same way as non-relativistic spin-0 (hypothetical) particles are described 
by the Schrodinger equation. The Dirac equation has been around for almost 80 
years and remains a formidable challenge to solve in all but the most trivial cases. 
Time-dependent problems are especially difficult. It is only with the use of modern 
computational resources and methods that situations of interest can be solved [4, 5]. 
These problems are very challenging, both theoretically and computationally. Their 
interpretation is also challenging, since they describe multi-particle situations, e.g., 
particle- ant iparticle pair creation. 
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1.1.1 Supercritical pair creation 

One of the predictions of the relativistic theory of the electron is that for fields of 
sufficient strength, so-called supercritical fields, the ground state will have its energy 
inside the negative-energy continuum. This is due to a feature of relativistic theories 
of particles: a lower continuum. The spectrum for a non-relativistic particle that 
vanishes asymptotically, such as the attractive Coulomb interaction contains both 
continuum states and bound states. The bound states start with the most negative 
energy (the rest energy is not included), namely the ground state and become denser 
as the energy rises (Rydberg series). At the ionization energy (the energy which 
the particle needs in order to become unbound) and above, there is a continuum of 
states the particle can occupy. These states, in the case of hydrogen, represent free 
electrons scattering from a proton. In Dirac theory this structure is also present, but 
there also exists a lower continuum below the bound states, starting at E = — mgC^ 
and extending to E ^ — oo. As the attractive potential becomes stronger, the 
system's ground state energy is lower, which causes no dramatic change in the non- 
relativistic theories. In Dirac theory, as the potential increases the ground state 
is lowered to eventually reside in the lower continuum (which is usually called the 
negative-energy continuum). This causes the ground state to change from a bound 
state to a resonant state, called a supercritical resonance state. 
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Resonances are formed, in general, when a bound state is embedded into a 
continuum of scattering states. Often this is due to the mixture of two potentials, 
e.g., when the Coulomb potential of an atom is perturbed by a unform electric field. 
In one direction, this results in the bound states having their eigenenergies among 
the continuum of scattering states. Thus, resonances are created (cf. section 2.3.1 
for more details). 

Supercritical resonances have all of the features of a usual atomic resonance with 
the primary difference being in the interpretation. The supercritical resonance is a 
pair creation resonance and not just the familiar scattering resonance from quantum 
mechanics. Supercritical resonances are characterized by two parameters in the 
usual way: the energy position and the lifetime. Since the lifetime is inversely 
proportional to the energy width, these parameters completely describe the shape 
of the resonance as it appears in the energy density of states (obtained from the 
spectrum) or in the scattering problem (viewed as scattering of the negative-energy 
continuum states). The resonance's shape is called a Breit-Wigner distribution 
which has the form 



where Ej-^s and T are the energy position and width respectively. It is, therefore, 
useful to determine the resonance parameters in order to facilitate this experimental 
search. 
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Two direct methods of obtaining the resonance parameters using the mapped 
Fourier grid method (cf. appendix A for details of the mapped Fourier grid method) 
are demonstrated in chapter 3. These extended methods are general enough to be 
applicable to other resonance phenomena in quantum mechanics. 

1.2 Analytic continuation methods to determine resonance 
parameters 

There are a few well-established methods for obtaining the resonance parameters 
using analytic continuation of the Hamiltonian (which means making the Hamilto- 
nian complex by an artificial transformation). These methods have been employed, 
compared, extended and improved in this thesis by using the supercritical reso- 
nance as a challenging example. This has led to very accurate determinations of 
the supercritical resonance parameters and shows general trends in the behavior of 
supercritical resonances. The improvements in the methods are equally applicable 
to normal scattering resonances in quantum mechanics. The resulting increase in 
performance is demonstrated along with the fact that several competing methods 
are shown to agree to very high precision. 

The method of complex scaling (CS) [6] and its extension, smooth exterior scal- 
ing (SES) [7, 8], are compared to the method of adding a polynomial complex 
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absorbing potential (CAP) [9, 10, 11], by determining the supercritical Dirac res- 
onance parameters. In this thesis, the CAP method is extended to the relativistic 
Dirac equation for the first time. All these methods use a parameter to charac- 
terize the transformation of the Hamiltonian from real to complex and then vary 
that parameter. Stability in the resonance parameters is then used to find the best 
approximation to the true resonance parameters; hence they are called stabilization 
methods. A comparison of several stabilization methods is carried out in chapter 4. 

A further extension of these analytic continuation methods was motivated by 
R. Lefebvre et al. in Ref. [12] where a Pade extrapolation was employed to increase 
the accuracy of the CAP method. In chapter 5, this idea is extended to SES. Con- 
vergent results are obtained by both methods when extrapolation is applied. The 
stability and parameter independence of the results are demonstrated. The extrap- 
olation also increases the accuracy of the analytic continuation methods, allowing 
for the exploration of effects beyond the first-order approximation to the two-center 
potential, which is used to model the potential of two highly charged ions. These 
investigations provide detailed knowledge of the behavior of supercritical resonance 
states, including how the mapped Fourier grid method handles supercritical res- 
onance states. Understanding the behavior of the mapped Fourier grid method 
with a supercritical resonance state is needed for time- dependent collisions with 
supercritical intermediate states. 
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In this thesis, the mapped Fourier grid method, a pseudospectral method, was 
chosen due to its suitabihty for the non-perturbative strong-field problems of in- 
terest. By extending the mapped Fourier grid method to the Dirac equation (cf. 
appendix A and Ref. [13]) we then have a suitable and efficient way to solve the 
Dirac equation numerically. The method builds a matrix representation of the 
Hamiltonian and obtains the eigenvalues and eigenvectors by diagonalization. The 
representation is carried out in coordinate-space and the method maps the coordi- 
nates allowing for an efficient coverage of phase-space [14]. It has been demonstrated 
to yield very accurate results for the hydrogenic problem, which is used as a test 
case. It is well suited for solving time-dependent problems due to the fact that the 
diagonalization yields a complete set of states. Expanding into this set of states, 
therefore, preserves unitarity. A key feature of the method, which is exploited in 
time-dependent problems (cf. chapter 6), is that the complete set of states contains 
the supercritical resonance. Thus, the resonance does not need to be added or 
constructed, as was done in prior calculations in the literature [1, 3]. 

1.3 Searching for supercritical resonances in collisions 

The collision of two heavy, fully-ionized atoms becomes supercritical when the 

atoms are a few tens of femtometers (10~^^m) apart. This is the only known 

route to supercritical resonance states that is accessible experimentally, making it 
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the focus of work seeking to demonstrate their existence. By treating each of the 
two nuclei in a collision as a homogeneously charged sphere separated by a distance 
R{t), a suitable model of the time-dependent electrostatic potential in the collision 
is obtained. The collision allows access to a supercritical potential, and therefore 
supercritical resonance states, by solving the time-dependent Dirac equation with 
this potential. Magnetic and retardation effects can be ignored, due to the slow- 
down of the nuclear motion near the Coulomb barrier. 

To solve the time-dependent Dirac equation, the propagator, which operates on 
an initial state and advances it according to the time-dependent hamiltonian, is 
approximated. The approximation assumes the system will not change over a short 
time interval. At, allowing for the use of the propagator with time-independent 
states. By solving the time-independent Dirac equation at a particular time, 
amounting to a fixed internuclear separation, R{t), and then using the approxi- 
mate propagator to move forward by At, the time-dependent Dirac equation for 
the collision is solved as a series of quasi-stationary steps. 

To propagate the state forward the eigenvalues and associated eigenvectors of 
the time-independent hamiltonian are needed. Also, in order to preserve unitarity, 
a complete set of states is required. As shown in Ref. [13], the mapped Fourier 
grid method can satisfy these conditions efficiently. Thus, by using the mapped 
Fourier grid method to solve the Dirac eigenvalue problem at successive internuclear 
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separations (which amounts to subsequent times) and by using these eigenstates in 
the approximate propagator to advance from t to t + At, it is possible to obtain the 
time evolution, not only for a single initial state, but for an entire set of independent 
initial states. 

1.4 Previous searches 

To demonstrate the existence of supercritical resonances, some signature of them 
must be sought. This work began in the 1970's and has continued to this day owing 
to the difficulty of both the theoretical predictions and the required experimental 
conditions. 

In 1981 the Frankfurt group published calculations of the resonance parameters 
and dynamical collision calculations. The group looked for evidence of supercrit- 
ical resonance states in the available data from the Gesellschaft fiir Schwerionen- 
forschung (GSI) [1]. The supercritical resonance calculations were performed in 
the center-of-mass frame for extended nuclei using the monopole approximation to 
the two-center potential. The group obtained results for the supercritical resonance 
parameters (to within a percent accuracy), and constructed the supercritical reso- 
nance state wavefunction for use in the collision calculation. This was achieved by 
solving the Dirac equation for a potential, which matched the two-center potential 

up to some small value, Vc, and then remained constant afterwards, with a value of 
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V{r) = V{rc). The choice of Vc was guided by the knowledge that the transition 
from bound-type to continuum-type character for the supercritical resonance state, 
took place around E^cs — V{r) ±mcC^ = (cf. section 2.3). This resulted in a wave- 
function that could be used to calculate the required matrix elements. Subsequently 
the resonance state was used together with an adiabatic quasi-molecular basis, in 
order to solve the coupled differential equations for the expansion coefficients and 
calculate the expected positron spectrum (and total positron creation probability). 
They found that although total positron production increases dramatically with the 
combined nuclear charge of different systems, the increase is smooth, making de- 
tection of the effects of the supercritical resonance decay difficult. The calculations 
were performed for different Fermi levels (the lowest level filled with electrons) to 
represent a collision in which a single /T— shell vacancy was created prior to closest 
approach. 

The GSI accelerator at the time was not able to provide fully stripped heavy 
nuclei and a solid, neutral, target was used to provide the collision partners. As 
expected, this damped the signal, since the electrons would be Pauli-blocked from 
being created below the Fermi level. 

Figure. 1.1 shows the results obtained for both supercritical collisions and col- 
lisions where even at closest approach the potential is not supercritical, i.e. a 
subcritical colhsion. 
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Figure 1.1: Positron spectrum from Reinhardt et al. [1] for collision systems 
with states up to SScx and 4Pi/20" occupied at a center-of-mass energy of 700MeV 
using a closed-coupling calculation. The results show no qualitative difference for 
supercritical collisions (U-U, U-Cf) as compared to subcritical collisions (Pb-U and 
Pb-Pb systems). 



Each curve represents the calculated positron spectrum for different collision 
systems from uranium-californium to lead-lead with states up to 3Sa and AVi/2<J 
filled initially. The lead-lead collision does not have a supercritical potential even at 
closest approach (^united ~ 164), but has the same qualitative shape as the uranium- 
californium spectrum (^united ~ 190) indicating that the decay of the supercritical 
resonance is not resulting in marked features. The authors compared their results to 
some unpublished data from GSI obtaining excellent agreement except for the case 
of uranium-uranium (experimental uranium-californium data were not available). 
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In 1983, the EPOS collaboration at GSI published their findings in Ref. [2] 
on a sharp peak (<80keV) from uranium-californium collisions at 6.05MeV/u. A 
similar peak was found in uranium-uranium collisions in unpublished work, but 
there were still questions as to whether the peaks had a common origin. Figure 1.2 
shows the collaboration's data and the peak in (a) is reported with a confidence 
level of <0.1%. The lines with error bars on both plots represent the positron 




200 /.OO EDO 800 1000 120D 



E,, [ keV ] 

Figure 1.2: Experimental positron spectrum from Schweppe et al. [2] for a 
uranium- californium collision system at 6.0MeV/u. The top plot is for prefer- 
entially backward elastic scattering 100° < 9 cm. < 130°, while the lower plot is for 
favored forward elastic scattering 50° < dc.m. < 80° for the ion (uranium). The 
dashed lines show the theoretical expectation for dynamical positrons including the 
nuclear background. Of interest is the peak in the top plot which is narrower than 
80keV. 
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spectrum, while the dashed hnes show the theoretically expected spectrum from 
dynamical positrons and the nuclear background. The top plot is for preferentially 
backward elastic scattering 100° < 9c,m. < 130° of the uranium ion, while on the 
lower plot it is for favored forward elastic scattering 50° < dc.m. < 80°. The authors 
initially thought this peak was from the spontaneous decay of the supercritical 
resonance, since it was near the expected energy of the supercritical resonance at 
closest approach. They argued that the possibility that the positrons were from a 
single internal pair conversion of nuclear transitions did not seem likely, since no 
emission line could account for the intensity of the peak. Explaining the peak as the 
decay of a supercritical resonance was a problem because to obtain such a narrow 
peak the nuclei would have to stick for a long time (> 10~^*^s). In another analysis 
of a very similar narrow peak in uranium-thorium collisions [15], the hypothesis of 
very long nuclear sticking, was found to be very unlikely. Instead, led by the fact 
that a similar peak was found in thorium-thorium and thorium-curium collisions, 
the authors examined the earlier hypothesis [16] that the peak was due to the decay 
of a previously unknown neutral particle. 

The theory group revisited the problem in Ref. [3] with the added detailed effects 
of nuclear sticking in an attempt to explain the narrow peak which showed up in 
the results at the EPOS and ORANGE experimental groups at GSI. The cause 
of the peaks are still under investigation [17] but it is now clear that the peak is 
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not due to supercritical resonance decay, since they were also found in subcritical 
collisions such as thorium-tantalum [18]. Miiller et al. were able to calculate the 
expected positron spectra for collisions in which the nuclei remain almost static 
at closest approach for up to 10 zepto-seconds (10~^'^s). Nuclear sticking times 
of a few zepto-seconds are experimentally obtainable in deep inelastic collisions as 
demonstrated in Ref. [19] and Ref. [20]. The sticking has a strong influence on 
the positron spectra, and it was shown that the peak of the positron spectrum 
will center closer to the resonance peak providing a clearer signal of the decay of a 
supercritical resonance shown in Fig. 1.3. 

These results also demonstrate that the sticking causes interference between the 
dynamical pair production, due to the motion of the highly charged nuclei, and the 
decay of the resonance. This is clearly visible from the second interference peak at 
around SOOkeV for the largest peak in the figure, which corresponds to a sticking 
time of 10 zepto-seconds. Calculations corresponding to a sticking time of 5 zepto- 
seconds (2°'^ curve from top) also exhibits a second peak although it is at around 
llOOkeV. The other dashed curves correspond to nuclear sticking times of 2 and 1 
zepto-seconds and still show a small shifting of each peak towards the resonance 
energy. This prior work remains the basis for comparison to the work presented 
here in chapter 6, and is the most recent (prior) work on the subject known to the 
author. 
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Figure 1.3: Positron spectrum from Miiller et al. [3] for bare uranium-uranium 
collisions at a center-of-mass energy of 740MeV. Dashed curves are for sticking 
times of 10,5,2,1 zepto-seconds in order from the tallest to the smallest peak. The 
solid line is for the collision without sticking. 

Work continues on these collisions, but it is mainly focused on explaining the 
sharp peaks and examining the possibility of longer nuclear sticking times; neverthe- 
less the search for supercritical states continues. The new GSI-FAIR experiments 
will perform bare uranium-uranium collisions removing the Pauli-blo eking of the 
previous experiments. The head-on collisions in these experiments represent the 
best hope for detecting supercritical resonances in the near future. 
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2 Theoretical background 



Using his equation, P.A.M. Dirac correctly predicted the existence of the positron 
more than 70 years ago and the equation has remained a focus of study ever since. 
The equation was intended to describe electrons relativistically, combining quantum 
mechanics with special relativity, but it ended up having a much broader scope. 



antiparticles in accord with relativistic kinematics. As a field equation, the Dirac 



forms the cornerstone of quantum electrodynamics (QED). QED is the most suc- 
cessful physical theory accounting for phenomena such as the electron's anomalous 
magnetic moment and the energy level shifts of hydrogen (Lamb shift) . 

Solutions to the Dirac equation are four-component wavef unctions, called spinors 
(or more precisely bispinors or Dirac spinors). In canonical Dirac-Pauli represen- 
tation, the Dirac equation for an electron (or positron) with mass is 



As a wave equation, the Dirac equation describes spin-| particles 




equation describes spin-i 



I fields. When coupled to quantized electromagnetism, it 




(2.1) 
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where the a matrices are commonly defined as 




\ 



(2.2) 



Here cxj are the 2x2 Pauh matrices and (3 = diag{l,l, —1, —1 



} with all other 



entries being 0. The solutions to Eq. 2.1 contain both positive- and negative-energy 



negative-energy solutions are a core feature of the Dirac equation. They hint that 
the interpretation of the Dirac equation as a single-particle theory is problematic. 
The interpretational difficulty posed by the negative-energy states is explored in 
the next section using hole theory. 

2.1 Hole theory 

The interpretational difficulties caused by the negative-energy solutions of the Dirac 
equation are due to the expectation that the relativistic theory of the electron 
should resemble the non-relativistic theory (the Schrodinger equation) and yield a 
probabilistic single-particle interpretation. The relativistic theory differs since it is 
valid at high energies allowing particles to be created and destroyed. Therefore, any 
description of the relativistic electron containing these phenomena will be a multi- 
particle theory. Figure. 2.1 illustrates the spectrum of the potential-free Dirac 
equation. The negative-energy states form a spectrum similar to the positive- 



four- component wavefunctions with both spin states ( 



I and — |). The existence of 
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energy continuum, but covering down to E 
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Figure 2.1: Illustration of the free Dirac equation spectrum. The upper continuum 
is similar to what is obtained in the free Schrodinger equation (except that the 
rest energy appears explicitly) but the lower continuum is only found in relativistic 
theories. The lower continuum is also called the Dirac sea and in hole theory is 
assumed to be completely filled with electrons preventing the decay of positive- 
energy particles into these states. The gap between the two continua is 2meC^, the 
minimum energy required to create a free electron-positron pair. 

pose the following problem: Positive-energy states could decay to negative-energy 
states continuously, thus emitting an arbitrary amount of energy [21]. This is 
not what is observed. Dirac proposed that the negative-energy states were all 
filled with electrons (the highest filled energy level is called the Fermi level). The 
Pauli principle then ensures that no positive-energy state could decay into a (filled) 
negative-energy state. This restored stability against radiative decay at the expense 
of losing the single-particle interpretation. Now all solutions to the Dirac equation 
contain information about the filled negative-energy states (the Dirac sea). They 
are inherently multi-particle solutions. In the low-energy regime the multi-particle 
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component is small and justifiably ignored. 



The Dirac sea allows the theory to incorporate the phenomenon of particle 
creation within so-called hole-theory. It takes into account the creation of particle- 
hole pairs (which can also be interpreted as particle-antiparticle pairs). In Fig. 2.2, a 
negative-energy electron is shown being excited into a state of positive energy. This 
can be interpreted as the creation of an electron and a positively charged particle. 
The absence of a negative charge in the Dirac sea will appear as a positively charged 
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Figure 2.2: Illustration of pair creation in hole-theory. A negative-energy electron 
receives sufficient energy to jump to a positive-energy state. The electron is then 
free and leaves behind a hole in the Dirac sea. This hole represents a missing 
negative charge and acts as a positively charged particle. The hole is the antiparticle 
of the electron since the electron can then fall back into the hole and fill it leaving 
the spectrum with no positive-energy electrons. 



particle, a hole. The excitation process of a negative-energy electron (c^^q) to a 
positive-energy state by a photon (7) is written schematically [22] as 
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(2.3) 



and is reinterpreted in hole-theory as 

7 ^ e;^>o + 4>o- (2-4) 

The reverse process, annihilation, where a positive-energy electron falls into an 
empty negative-energy state filling a hole (vacancy) is written schematically as 

e£>o ^ ^E<o + 27, (2.5) 
which is reinterpreted in hole-theory as 

e£>o + 4>o 27. (2.6) 

The positively charged particle with positive-energy, e^^Q, is therefore to be inter- 
preted as a positron since it annihilates the electron. To do this the hole must be 
viewed as a particle, which means that its dynamical quantities are those of a hole 
with positive-energy. For instance, the spin of the positron must be the opposite of 
the missing negative-energy electron, since a missing spin-up electron would behave 
as a spin-down positron and vice versa. The momentum must also change sign. The 
absence of momentum p in the Dirac sea will appear as the presence of a particle 
of momentum —p. In this way, hole theory allows for the Dirac equation to de- 
scribe physical situations in which the number of particles is not conserved making 
it much more useful than a single-particle theory. Dirac predicted the existence of 
the positron using his equation, despite initially not being taken too seriously as 

demonstrated by this quote from W. Pauli's Handhuch article in Ref. [22]: 
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Recently Dirac attempted the explanation, already discussed by 
Oppenheimer, of identifying the holes with anti-electrons, particles of 
charge +|e| and the electron mass. Likewise, in addition to protons 
there must be antiprotons. The experimental absence of such particles 
is then traced back to a special initial state in which only one of the 
two kinds of particles is present. We see that this already appears to be 
unsatisfactory because the laws of nature in this theory with respect to 
electrons and antielectrons are exactly symmetrical. Thus 7-ray photons 
(at least two in order to satisfy the laws of conservation of energy and 
momentum) must be able to transform, by themselves, into an electron 
and an antielectron. We do not believe, therefore, that this explanation 
can be seriously considered. 

This article ended up appearing in print after C. D. Anderson had already verified 
the existence of the positron after analyzing particle tracks and looking for a particle 
with the electron's mass and opposite charge [23]. 



2.2 The Dirac equation with a potential 

To discuss more than free particles a potential should be added to the Dirac equation 
and there are some interesting consequences with respect to the negative-energy 
states. The relativistic covariance of the Dirac equation has the feature that a 
scalar potential will affect the negative- and positive-energy states identically. For 
details of the behavior of the Dirac equation with a scalar potential the reader is 
referred to Ref. [21]. 

A vector potential of relevance to the present work, is the electrostatic Coulomb 
potential. The Coulomb interaction is the time-component of the four-vector po- 
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tential in the Coulomb gauge. It is written in natural units {h = = c = 1) 
as 



where Z is the charge and a is the fine-structure constant. The Dirac equation with 
an attractive Coulomb potential, the relativistic hydrogen atom, represents a major 
success in the spectroscopy of atomic hydrogen. It accounts for the corrections due 
to L ■ S coupling, i.e., orbital angular momentum coupling to the intrinsic spin an- 
gular momentum of the electron. This is due to the spin being a fundamental piece 
of the Dirac equation. As an aside, we mention that high-precision spectroscopy 
revealed an inadequacy of the equation as early as 1947: the Dirac equation pre- 
dicts degeneracies in energy levels (e.g. 25*1 vs 2Pi) which are only approximate 
in nature. The 2Pi state can decay radiatively to the ground state by one-photon 
emission, and thus becomes a resonance in quantum electrodynamics. In the Dirac 
theory all excited states are infinity long lived (exact energy eigenstates). 

A key difference between scalar and vector potentials is that a vector potential 
acts differently for positive- and negative-energy states. The relativistic hydrogen 
atom has the same continuum spectrum as the potential-free Dirac equation with 
added states in the gap just below the positive continuum shown in Fig. 2.3. These 
states are called bound states since they are localized in coordinate-space. Also de- 
picted in Fig. 2.3 is the spectrum of hydrogen-like uranium (U^^^), although not to 
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Figure 2.3: Illustration of the hydrogen spectrum for hydrogen (left) and hydrogen- 
like uranium (right) with neither to scale. The lowest bound state, the ground state, 
is lower for hydrogen-like uranium than for hydrogen. Thus it takes more energy 
to ionize the ground state electron bound to a uranium nucleus than to a proton. 

scale. The state closest to the negative-energy continuum is the ground state since 
it is the lowest positive-energy state. It is called the IS state in atomic systems 
(where the digit stands for the principle quantum number ranging from 1 to infin- 
ity and the letter stands for the orbital angular momentum with S representing no 
orbital angular momentum or a spherically symmetric state [24]). As the strength 
of the central binding potential (such as the Coulomb potential of a nucleus) be- 
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comes stronger the ground state, as well as the other bound states, will sink deeper 
into the gap. This explains why it takes more energy to excite the ground-state 
electron to the continuum in hydrogen-like uranium (compared to hydrogen). In 
the Schrodinger theory this energy scales like Z^. 

Solving the Dirac equation for the Coulomb potential, as found in Ref. [22] or 
[21], one obtains the Sommerfeld formula for the energy of the bound states. 



E= , (2.8) 



where Z is the charge of the nucleus, is the relativistic principal quantum number 
ranging from to infinity (where the principle quantum numbers are related by 
n = rir + {j + |)), a is the fine-structure constant and j is the total angular 
momentum quantum number and is equal to | for S-states. As Z increases, the 
denominator in Eq. 2.8 increases. A state with the same rir and j will therefore 
have a lower energy for higher Z. A lower energy eigenvalue means the electron will 
be more tightly bound to the nucleus. Figure 2.3 shows this difference for Z = 1 
(hydrogen) and Z = 92 (hydrogen-like uranium). 

While uranium-238 is already one of the heaviest stable nuclei, the following 
question may be posed: What are the consequences of having a (hypothetical) 
nucleus with higher charge Z7 A pure Coulomb potential originates from a point 
source and the Dirac equation for this case does not yield solutions for a charge 
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larger than Z > 137 (for Z < 137 one obtains a ground-state energy very close to 
the center of the gap, i.e. Eis ~ 0). 

The pure Coulomb potential is unrealistic for heavy nuclei. A more accurate 
model of a nucleus involves a finite charge distribution (such as a homogeneously 
charged sphere) which allows solutions for any Z value. For very large Z (Z > 
170) the ground state energy resides in the negative-energy continuum leading to 
interesting effects discussed in the next section. Unfortunately, there is no evidence 
for the existence of such very heavy nucleus {Z > 170) in nature. The effects are 
mimicked in fully-ionized atomic collisions when two nuclei are nearly touching. 
We, therefore, will use the model of a single heavy nuclei as a simple explanatory 
tool, but look for the effects in the collisions of heavy, fully-ionized atoms. 

2.2.1 Supercritical ground states 

When considering the difference between the two spectra in Fig. 2.3 and examining 
the Sommerfeld formula (cf. Eq. 2.8), one can see that increasing the charge of the 
nucleus will result in deeper bound states. The ground state is always the most 
affected and decreases fastest in energy. Eventually, as the charge of the nucleus 
is increased beyond about 170 (depending on the model of the nucleus used), the 
ground-state energy will reside in the negative-energy continuum. An electron 
bound in such a state would require a photon with sufficient energy to create an 
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electron-positron pair [E > 2mcC^) to ionize the electron. If this state was empty it 
would be possible for one of the negative-energy electrons, with E ^ Eis, to fill the 
ground state. This is because the empty ground state introduces a hole (which is 
initially "bound" to the nucleus) into the negative-energy continuum, since unlike 
the negative-energy states, this ground state is still a partially localized state. Any 
electron which occupies it will remain very tightly bound since it is Pauli-blocked 
from tunneling out to the negative energy states. Once a negative-energy electron 
fills the ground state the hole will then be in the negative continuum. This hole 
will escape from the nucleus since it acts as a positively charged particle and is 
repelled by the highly charged positive nucleus. This situation can be reinterpreted 
as pair production, since the previously fully-ionized nucleus now has an electron 
in the ground state and a free positron has escaped (this is commonly referred to 
as bound-free pair production, since the electron is bound and the positron is free 
in the final state). The hole (positron) will have an energy very close to that of the 
ground state since it is most likely that the negative-energy electron which filled 
the ground state has the same energy. 

It is also possible for the potential to be even stronger causing the ground 
state to become more deeply embedded in the negative-energy continuum. This 
increased potential will also lower the energy of the other bound states. The 2Pi, 
the first excited state, becomes supercritical at Z > 180 [25]. This will result in 
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other states, besides the ground state, becoming supercritical These states will 
behave as already described and, therefore, although the description given is for a 
supercritical ground state they equally apply to any bound state which becomes 
supercritical. 

The spontaneous pair production due to the potential of such heavy nuclei is 
often referred to as charged- vacuum decay [25]. This is because the vacuum itself 
is unstable and decays by spontaneous pair production. This is reminiscent of the 
Klein-paradox (where a step potential of Vq > 2mcC^ shows reflection and transition 
of an impinging relativistic wavepacket) and a feature of the multi-particle nature 
of the Dirac equation with very strong potentials [5]. While such large potentials, so 
called supercritical potentials, are not available from single nuclei, colliding nuclei 
can come into close contact and for a short time create a potential large enough 
for charged-vacuum decay to take place. It may also soon become possible for an 
intense laser pulse to add to the potential of a nucleus which would also allow for 
charged-vacuum decay. 

2.3 Resonances and supercritical resonances 

When the energy of a bound state is in the negative-energy continuum, becoming 

a supercritical resonance state, its wavefunction will change its character. A bound 

state is characterized by a wavefunction which is localized around the nucleus and 

30 



decays exponentially at larger distances. A continuum state is non-localized and 
always oscillating. Both types of states are shown in Fig. 2.4 with the bound state 
being the ground state of uranium and the continuum state E = l.SmeC^. When 
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Figure 2.4: The left plot is the radial probability density (|r\l'(r)p) of the uranium 
ground state, ISi, as function of the radial coordinate (units h/(mec) ~ 386fm). 
The right plot is the radial probability density of continuum state with E = l.SmeC^. 
The bound state is localized since the probability density is concentrated near the 
origin (the nucleus) and decays exponentially with increasing radial distance. The 
continuum state oscillates for all distances since it characterizes a non-local state. 

a state becomes supercritical it takes on the characteristics of both a bound and a 
continuum state and is called a supercritical resonance. 

2.3.1 Resonances 

A resonance state in the current context is defined as: a scattering system which 
has sufficient energy to break up into two or more subsystems, if it is long-lived 
compared to the collision time [26, 27]. Illustrated in Fig. 2.5 are the three possibil- 
ities of the scattering of two oppositely charged bodies. The first situation results 

when the projectile does not have sufficient energy to escape and becomes bound 
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Figure 2.5: The three possibihties of scattering two oppositely charged particles. 
The left hand illustration represents the projectile becoming bound. The middle 
illustration shows the projectile being deflected and the third is a mixture of the 
two and called a resonant state. 

to the target. The second possibility occurs when the projectile is simply scattered 
by the target, and the third is a mixture of the two. In the last case the projectile 
becomes bound for a time which is long compared to the collision time and then 
scatters away. The lifetime of this quasi-bound state is defined by the amount 
of time the system will stay bound before having a 63% probability (1 — e~^) of 
decaying into two separate parts (projectile and target). 

When scattering electrons with different energies from a target with a resonance 
at some energy, E^-es, a peak will show up in the scattering cross section, cr{E) at 
E ^ i^rcs- This peak will have a width, F, which is inversely proportional to the 
lifetime, r, of the resonance by r = h/T. It will take the shape of a Breit-Wigner 
distribution (cf. Eq. 1.1). The narrower the peak the longer the projectile electron 
will stay bound to the target before scattering away. A common example in atomic 
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physics, used in explaining high-harmonic generation, is the resonance formed when 
the Coulomb potential of an atom is perturbed by a strong electric field (from a 
laser for instance), illustrated in Fig. 2.6. Electrons scattered off an atom with the 

E=0 E=Eo 

Figure 2.6: Illustration of the radial Coulomb potential of an unperturbed atom 
(left) and an atom in the presence of a strong external electric field (right). The 
potential (along some axis) is the black line and the nucleus is depicted as a solid 
sphere. The ground state wavefunction is shown in blue with it being a bound state 
when the electric field is off and a resonance state when the field is on. 

external electric field will resonate if their energy is of the order r/2 away from 
Ej-cs or less. This can be thought of as the projectile electron tunneling through the 
potential barrier to the bound region, remaining in the quasi-bound state for a time 
and then tunneling back out. Resonance states are therefore different from both a 
bound or regular continuum state and their time-independent wavefunctions reflect 
this by being a mixture of the two. A resonance state wavefunction has a bound 
part which is localized around the potential well (deepest part of the potential), but 
instead of exponentially decaying the tail oscillates as a continuum state as shown 
in the right side of the illustration in Fig. 2.6. The tail gives some amplitude to 
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the resonance state away from the nucleus where the continuum states have their 
amphtude. The matrix elements (which determine the connection between different 
states by (\l/i(r)|/(r)|\l/2(^)) = / ^i('")/('")^2(^) for some coupling /(r) such as 
an external static electric field) between the resonance state and the continuum 
states will thus be larger than continuum-bound matrix elements. In this way, the 
continuum (or scattering) states are connected to the bound part of the resonance 
state. 

2.3.2 Supercritical resonances 

Supercritical resonance states have the same form for their wavefunctions, but the 
interpretation differs. In the case of a single supercritical nucleus a ground-state 
vacancy can decay into a bound electron and a free positron. This takes place 
despite there being no externally scattered particles. The situation may instead 
be viewed as the scattering of the negative-energy electrons from the supercritical 
vacancy state. The behavior is then similar to electrons scattering off an atom in an 
external electric field. The supercritical resonance has a lifetime and the resonance 
peak will be at the energy of the supercritical bound state. The interpretation is 
different, since in the case of supercritical resonances the outgoing state will be a 
bound-free pair from an initial vacuum state. 

The change from bound state to supercritical resonance state results in the 
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exponentially decaying tail being replaced with an oscillating one. The lifetime of 
the vacancy thus becomes finite. When the supercritical resonance state is occupied 
it is stable. This is due to Pauli-blocking which prevents any tunneling to the (filled) 
negative-energy continuum. As expected, this atom is stable just as a subcritical 
atom (an atom which has no supercritical states). The lifetime refers only to an 
unoccupied supercritical state which then decays by pair creation. This is analogous 
to the atom in an external electric field since if the resonance state is filled other 
electrons cannot tunnel in (until the resonance electron tunnels out). 

2.3.3 Cause of the supercritical resonance 

An atom perturbed by an external electric field creates a resonance state by modi- 
fication of the pure Coulomb potential. This creates a finite-size barrier, depicted 
in Fig. 2.6. Three regions result: a bound region near the atom, a region where the 
potential barrier is high, and an external region where the potential is weak. The 
region where the potential is high but of finite extent is sometimes called the barrier 
region, since this is where the electron can tunnel into or out of the resonance state. 
It is also the region where the wavefunction crosses the border from being a bound 
state to an oscillating continuum-type state. In most situations in atomic physics 
where resonances occur it is due to a potential barrier being formed by the mixture 
of two or more potentials. 
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In the case of supercritical resonances no other potential is present except for 
the modified Coulomb potential (called modified due to the use of a finite-size 
nucleus). It is not the finite nucleus which is the cause of the potential barrier, since 
the transition from within the nucleus to outside the nucleus must be continuous. 
Instead, it is, in fact, the gap between the positive and negative continuum states 
which acts as the tunneling region and is illustrated in Fig. 2.7 for a subcritical and 
supercritical ground state. 



Figure 2.7: In black the continuum limits, given by V{r)±mcc'^, are shown as a 
function of the radial distance. The probability density, |r\l/(r)p, of the ground 
state for each situation is sketched in blue and displaced vertically using its energy 
as a baseline. The subcritical ground state has its tail in the gap and thus decays 
exponentially. The supercritical ground state is more deeply bound and its tail 
resides in the negative-energy continuum allowing it to oscillate. 

The illustration shows the ground state vertically displaced by its energy and 
the energy spectrum is given as a function of the radial distance perturbed by the 
Coulomb potential. Near the origin the spectrum is highly perturbed but further 




away less so. In the case of a subcritical potential the tail of the wavefunction is 
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completely contained within the gap between the positive and negative continua. 
When the potential is supercritical the spectra are further perturbed causing the 
binding energy of the ground state to increase. This results in the ground state's 
lobe compressing towards the origin. It also decreases the energy of the state, 
lowering it in the spectrum. Consequently the tail of the ground state ends up 
inside the negative-energy continuum. The tail of the ground state can therefore 
oscillate since it is not in the gap. 

The ground state now has the characteristics of a resonance state. The resonance 
state is, therefore, the result of having two continua separated by a gap and a vector 
potential, not by the interplay of two separate potentials. 

2.4 Accessing supercritical resonances 

Stable nuclei with sufficient charge to have a supercritical ground state are unlikely 

to exist. If they could be created experimentally, their decay time would very likely 

be much shorter than the decay time of the supercritical ground state which is found 

to be on the order of 100 zepto-seconds (100xlO~^^s). Supercritical resonances 

are predicted by the Dirac equation making them a fundamental prediction of 

the relativistic theory of the electron. Therefore, other avenues must be sought 

to verify this prediction. Currently, the most promising approach is to use fully 

ionized atomic collisions (see section 1.4 for previous attempts). 
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Atomic collisions with fully ionized heavy atoms yield a supercritical ground 
state when the nuclei are very close, since the total potential can become supercrit- 
ical. To visualize the situation it is helpful to use the quasi-molecular picture where 
the two nuclei are considered to form a slowly varying quasi-molecule with changing 
internuclear separation. The system then has a single (molecular) spectrum. The 
potential of the two nuclei is that of two displaced and moving modified Coulomb 
potentials. At any instant in time, the ground state is then a quasi-molecular 
ground state (called the ISo" where the additional a states the radial wavefunction 
is symmetric) which covers both nuclei. While the potential is no longer spherically 
symmetric or time-independent, it becomes experimentally accessible and yields 
supercritical resonance states that may be detectable. 

2.4.1 Dynamical and resonant pair production 

In fully-ionized collisions with a supercritical resonance, two types of pair creation 
will be present. The first type is due to the spontaneous decay of the supercritical 
resonance, as has been described in the previous section. These electron-positron 
pairs are not directly due to the dynamics of the collision. The free positrons (and 
ISo" electrons) are created due to the instability of the vacuum in the presence of 
an intense field; within Dirac hole theory the process is represented by the decay 
of the supercritical resonance state. The process is, therefore, called resonant or 
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spontaneous pair creation. The other type of electron-positron pair creation is what 

is normally observed in heavy-ion collisions. It is due to the changing potential 

caused by the moving nuclei. It is referred to as dynamical pair creation. 

Hole theory provides a useful picture of the situations as illustrated in Fig. 2.8. 
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Figure 2.8: Illustrative explanation of dynamical pair production (left) and resonant 
pair production (right). Only a few of the infinite continuum states are shown. 
The circles represent holes while the filled circles represent electrons. Dynamical 
pair production occurs when a negative-energy electron absorbs some of the kinetic 
energy from the motion of the nuclei and is excited to a positive-energy state leaving 
a vacancy behind. Resonance pair production occurs when the vacant ground state 
joins the negative-energy continuum which allows a negative-energy electron with 
E ^ Ej-cs to tunnel in and become bound to the nucleus while the vacancy, now in 
a negative-energy continuum state, escapes. 



The left picture shows dynamical pair creation. In hole theory, this amounts 



to a negative-energy electron absorbing some of the colliding nuclei's kinetic en- 
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ergy and jumping into a positive-energy state. The electron leaves a hole in the 
negative-energy continuum representing a free positron. The electron can jump to 
a bound state creating a bound-free pair or, with sufficient energy, to a continuum 
state representing free-free electron-positron pair production. The energy for this 
transition comes from the kinetic energy of the collision. Since heavy-ion collisions 
are generally performed at center-of-mass energies well above 200MeV, the l-2MeV 
required for this process is small , i.e., it is readily available. In collisions where par- 
tially ionized atoms are used the lower bound-free channels require the inner-shell 
electrons to be excited or else the channel is Pauli-blocked, severely suppressing the 
channel. 

As the nuclei approach each other, the quasi-molecular ground state of the sys- 
tem is lowered. A negative-energy electron would require less energy to be excited 
to a bound state. This is why bound-free pair production, or most importantly 
ground state-free pair production, is the dominant channel. In cases of collisions 
which involve supercritical resonances, dynamical pair creation is further enhanced. 
This is again due to the smaller amount of energy required to excite negative-energy 
states into the (now lower) bound states and into the supercritical ground state. 
Energetically it may be expected that the negative-energy states closest to the gap 
would be those most likely to be excited by the collision into a positive-energy 
state, but this is not the case. The continuum states near the gap are states of low 
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momentum. Thus, they cannot react to the colhsion or connect to the locahzed, 
low-lying, bound states during the collision. 

The second type of pair production, illustrated in Fig. 2.8, has been discussed 
in the previous section, the decay of a supercritical resonance. During the time in 
which the system has a supercritical ground state (or possibly other states as well) 
the bound hole can tunnel out into the negative-energy continuum, or equivalently, a 
negative-energy electron can tunnel into the bound region of the supercritical state. 
Clearly, negative-energy states of similar energy to the supercritical resonance are 
the ones which can tunnel into it, and the resulting positrons will be close (on the 
order of r/2 away) to the energy of the supercritical resonance. This type of pair 
creation is exclusively bound-free since the electron is bound to the nuclei once it 
tunnels in from the negative-energy states. 

These two types of pair creation, therefore, have significant overlap since the 
exclusive resonant pair creation adds to the dominant channel for dynamical pair 
creation, ground state-free pair production. These two effects interfere with each 
other making the determination of the single cause (dynamic or resonant) of the 
electron-positron pair production difficult to discern. An important difference is in 
the positron spectrum. While the electrons will primarily be in the ground state, 
since this is the dominant channel for dynamical pair creation and the only channel 
for resonant pair production, the positrons from resonant pair production will be 
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peaked around the energy of the supercritical resonance. This peak will be smeared, 
since unlike in the static case the supercritical ground state will change due to the 
nuclear motion. As the nuclei approach, the state's energy will sink deeper in 
the negative continuum and then come back up again as the nuclei separate. The 
resonant pair production will then be a smeared version of all the intermediate 
supercritical resonances. The smearing is not equal though, since the deeper the 
resonance the more momentum the negative-energy states have resulting in shorter 
decay times of the state. The shorter decay time implies that the resonant pair 
production will be dominated by the deepest supercritical resonance, created at 
the closest approach of the nuclei. 

The supercritical resonance occurring at closest approach for currently exper- 
imentally available fully-ionized atoms has a lifetime of the order of 10~^^s. The 
interaction time in these collisions at closest approach is of the order of 10~^^s, 
resulting in a very small signal for the resonance production peak. Maximizing the 
time the system is supercritical is, therefore, of the utmost importance, if the goal 
is to observe charged vacuum decay. 

2.4.2 The trajectory 

The main objective of this work is to test the prediction of the Dirac equation with 
regard to supercritical resonances; consequently we are interested in trajectories 

42 



which try to maximize the time in which the system is supercritical. To this end a 
Coulomb trajectory is chosen. 

A Coulomb trajectory for heavy ions involves relatively slow nuclear motion 
('^^max ~ 0.11c for bare uranium-uranium) compared with the velocity of a ground 
state electron. The center-of-mass energy is just enough for the nuclei to slow 
down as they get closer due to the mutual Coulomb repulsion, and come to a 
stop at closest approach before being accelerated away. The Coulomb trajectory is 
obtained by solving the classical force problem due to the repulsion of two nuclei 
of charges, Zi and Z2 given by 

,.f = (2^9) 

in natural units where /i is the reduced mass of the system. The Coulomb trajectory 
is illustrated in Fig. 2.9 where the arrows represent the velocity and show the ions 
slowing down as they approach each other. 
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Figure 2.9: An illustration of the Coulomb trajectory for two identical fully-ionized 
atoms from left to right (using six snapshots). The arrows beside the nuclei rep- 
resent the velocity showing how the nuclei slow down as they approach due to 
the Coulomb repulsion. The green arrow shows the distance of closest approach, 
determined by the collision energy. After reaching closest approach resonant pair 
creation is depicted with the green electron being bound in the ground state of one 
of the nuclei and the yellow positron scattering away. 

This trajectory offers the advantage that the nuclei are traveling very slowly 
when the system is supercritical, reducing the dynamical pair production and max- 
imizing the time in which the supercritical resonance has to decay. This differs 
from a straight-line trajectory where projectiles are of sufficient kinetic energy to 
have very little deviation due the collision. The Coulomb trajectory in contrast has 
a much lower center-of-mass energy reducing the amount of dynamical positron 
production. The ratio of the two pair production channels is therefore optimized 
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since the dynamical production acts effectively as a background to the supercritical 
signal. 

2.4.3 The reference frame 

The full Dirac equation is covariant yielding the same results regardless of the 
reference frame. However, the approximations which are introduced to solve the 
collision render the results reference-frame dependent. There are two convenient 
reference frames for the collision. 

In the target frame, one chooses one of the nuclei, often the heaviest, as the 
origin. The other nucleus is then treated as a projectile. The projectile's potential 
is then a perturbation to the modified Coulomb potential of the target. This intro- 
duces the time dependence into the system (due to the movement of the projectile). 
This reference frame has the advantage that the initial states can be taken as sim- 
ple atomic states since at t — > — oo the system is that of a single atom. The final 
states can also be projected onto the single atomic basis if propagation for t — > cxo 
is impractical. 

The center-of-mass frame on the other hand treats both nuclei equally, a desir- 
able feature when the nuclei are close and of similar charge. Since the eigenvalues 
of the states are frame dependent (the eigenvalue is the time component of a four- 
vector) the center-of-mass frame will have deeper bound states. The supercritical 

45 



resonance states will also be deeper resulting in larger positron yields. Therefore, 
care must be taken when choosing the reference frame to ensure meaningful results. 

2.5 Nuclear sticking 

Collisions where the nuclei end up within a nuclear radius of the each other result in 
deformations of the nuclei. This is due to the Coulomb repulsion. There is also the 
possibility of some effects due to the nuclear force. In deep inelastic heavy-ion col- 
lision experiments with lead-uranium, uranium-uranium and uranium-californium 
collisions, detailed trajectory reconstruction found that the nuclei may have re- 
mained in contact for around a zepto-second [19]. The causes of this effect are 
still unclear. Some groups predict that a sticking time of up to 10 zepto-seconds 
may be achievable [17, 28, 29]. This would give a strong enhancement to the su- 
percritical resonance decay signal due to the increased time for the supercritical 
resonance to decay (cf. Fig. 1.3). This would provide strong evidence for resonant 
pair production as well as long nuclear sticking times. 

2.6 Modern interpretation of hole theory 

Hole theory has been used in this section to understand the cause of the supercritical 
resonance and what it is, but a common question arises: what is the validity of hole 
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theory? In the current context hole theory provides a meaningful way to understand 
the negative-energy states and multi-particle aspects of the theory, but it is only 
an explanatory tool for the electron-positron field. In the modern view, it has 
been supplanted by quantum electrodynamics (QED) which deals properly with 
multiple particles and has only positive-energy states. An example of the failure of 
hole theory is simple to illustrate using a phenomenon beyond the electron-positron 
field. Any phenomenon in which the electron or positron numbers are not conserved 
would suffice. For instance P~^-dec&y, 

p^n + e^ + (2.10) 

where a proton, p, decays into a neutron, n, a positron, e"*", and a neutrino, v. 
Here a positron was created, a vacancy in hole theory, without an electron being 
created. This is expected since hole theory is not valid beyond the electromagnetic 
field interaction and the above process is a weak field decay. 

In the present collision work, the final quantities of interest, discussed in sec- 
tion 6.1.1, are derived from QED as in Ref. [30] and in appendix B for a non- 
orthogonal basis. The quantities necessitate that all initial Dirac wavefunctions be 
propagated through the collision in order to obtain the desired positron and elec- 
tron creation spectra in order to search for evidence of a supercritical resonance. 
Therefore, while hole theory provides some guidance to understanding the super- 
critical resonance state, it is ultimately QED which directs the search and provides 
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the quantities needed. 
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3 Time-independent supercritical states 



In single-particle (non-relativistic) quantum mechanics, bound states become reso- 
nances when the potential is modified such that the previously bound particle can 
escape to infinity (e.g., when the linear Stark potential is added to the Coulomb 
interaction). More generally, a system which has sufficient energy to break up into 
two or more subsystems is called resonant if it is long lived compared to the collision 
time [26, 27]. Such a state is described by the mean energy position i^ros (which 
is usually shifted from the bound-state eigenenergy) , and by the lifetime r. The 
latter provides a measure of the decay time of an initially prepared quasi-bound 
state which undergoes exponential decay. 

In scattering theory, the total cross section resonates (reaches the unitarity 
limit) if the particle energy is close to Ej-^s and can be described by a Breit-Wigner 
distribution (cf. Eq. 1.1) which is characterized by the energy width parameter 
r. According to the energy-time uncertainty principle, the temporal behavior of 
a decaying bound state and the cross section shape are related by Fr ~ /i. A 
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determination of T and Ej-^s and a fit to the Breit-Wigner shape function normally 
requires a solution of the scattering problem for many energies in the vicinity of 



Resonances embedded in the positive continuum are characterized by E-^ = 
E'rcs — Supercritical resonance energies however are given by E^ = Ej.cs + iT/2. 
This is understood by reinterpreting the supercritical resonance for an electron with 
negative energy as a positive-energy positron resonance propagating backwards in 
time by CPT symmetry. The time reversal transformation necessitates a positive 
imaginary part of the energy eigenvalue to ensure that the state decays as time 
propagates to negative infinity. This follows since the time-dependent part of the 
wavef unction is given by 



This time dependence then ensures that the total probability in the negative-energy 
hole state decays exponentially with time scale H/T as t ^ — oo. 

We note that some of the work in chapters 3 and 4 was originally published in 
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[31]. 
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3.1 A supercritical potential 



We use a physically relevant model of two extended nuclei separated by distances 
of the order of a few nuclear radii. The nuclear radius is taken to be = 1.2\^ 
fm, where A is the number of nucleons [1]. The nuclei are assumed to be displaced 
homogeneously charged spheres with a separation of R between the charge centers. 
In the center-of-mass frame the monopole potential for each of the two nuclei with 
respective charge Zj and radius i?^*^ is given (in units of^ = c = mc = l,Z = 
Z„/?n = /2«)by 
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where r± = R/2 ± R^. The expression is obtained by using the potential for a 
homogeneously charged sphere at the origin and displacing it by R along the z-axis 



and expanding it into Legendre polynomials, 
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(3.3) 



By inverting the sum in Eq. 3.3 using the orthogonality relation of Legendre poly- 
nomials, J^^ Pmix)Pn{x)dx = — 27^^*^"'™' "^'^ obtained for Vo(r) for each 
nuclei. 
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3.2 Supercritical states using the mapped Fourier grid method 

The mapped Fourier grid (MFG) method provides approximate resonance param- 
eters directly by two separate means. These methods also illuminate the behavior 
of the MFG method when a supercritical state is present, which is particularly im- 
portant for time- dependent heavy-ion collisions in which supercritical intermediate 
states strongly influence the time evolution of the matter field. 

3.2.1 The projection method 

Even though it is unlikely that a discretized continuum eigenstate of the Hamilto- 
nian matrix will fall exactly on the mean energy of the resonance, the closer such a 
state is to the mean energy the more it resembles the resonance state. The bound 
(small-distance) part of a supercritical IScr resonant state is similar to a subcritical 
ISo" state (although with a somewhat compressed lobe). The inner product of a 
normalized subcritical and supercritical IScr equals approximately one. Thus, the 
projection of a subcritical IScr onto the eigenbasis of a supercritical hamiltonian 
quantifies how similar the latter states are to the ground state. Therefore, by taking 
the inner product of each supercritical eigenstate with the subcritical ISa state one 
computes a probability density vs energy, E^, based upon \{(f){E,^)\ipisa)\'^- This 
probability distribution displays a Breit-Wigner shape. By fitting to the Breit- 
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Wigner shape (cf. Eq. 1.1), the parameters of the supercritical resonance {Ej.cs,T) 
are obtained. This method does not need to have any background subtracted since 
it is flat away from i^rcs- The result is, however, somewhat sensitive to the range of 
energies used in the fit and needs enough states to properly cover the distribution. 

3.2.2 Supercritical resonance states in the density of states 

A further method from which the supercritical resonance parameters can be esti- 
mated is based upon the density of states from a diagonalization of the discretized 
Dirac Hamiltonian. The density of states (in arbitrary units) 



is calculated from the quasi-continuum eigenvalues, E^. The fit requires the ap- 
proximation to the local background density, and, thus, the results are somewhat 
sensitive to the form of the fit. Using a single rational term for the background 
with a Breit- Wigner distribution added the density of states can be fitted by 




(3.4) 
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+ A2 
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(3.5) 




where Ai, A2 and q are fit constants. 
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3.2.3 Results obtained using the mapped Fourier grid method 

The success of the aforementioned methods is demonstrated for the example of 
the U^^^-Cf^^"*" quasi-molecular system for an internuclear separation of R = 20fm 
where Au = 238 and Acf = 251. The resuhs of 2 |(0(E^)|7/;is^(R = 45fm))|^ /{E^ - 
Ey^i) are shown as (red) crosses in Fig. 3.1 (the factor of IjiEy — E^^i) is due 
to the continuum states being quasi-continuum states and having a width). The 
(f>{E^) states were obtained by solving the Dirac hamiltonian for the potential in 
Eq. 3.2 with R = 20fm. Fitting the results in the range of the resonance to a 
Breit-Wigner the (blue) line in Fig. 3.1 is obtained with the resonance parameters: 
Ercs = -1.75795 meC^ and F = 4.12 keV. 

Building the density of states using Eq. 3.4 for the same supercritical system 
the (red) crosses in Fig. 3.2 are obtained. Fitting the results with Eq. 3.5 results 
in the (blue) line. The fit yields E^^s = — 1.75794 mgC^ and F = 4.17 keV for the 
resonance parameters. 

These results compare well to those of Ref. [1] which obtained values of E^^s = 
— 1.763meC^, F = 4.12keV using phase-shift analysis and Ej-^s = — 1.761meC^, 
F = 4.10keV using the truncated potential method (described in section 1.4). The 
projection method, while sensitive to the energy range used for the fitting, ob- 
tains a width closer to these reference values, but a slightly lower E^^s position. 
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Figure 3.1: The projections, 2 \{(j){E^)\ipisa)f /{E^ — Ey_i), of a bound ^/^iso- state 
onto the supercritical U-Cf quasi-continuum states (p{E^) (calculated for a = 
3900 and s = 6000 basis) are shown as crosses (+). The subcritical bound ipisa was 
obtained for a separation of i? = 45 fm, while the supercritical basis is for i? = 20 
fm. The line represents a Breit-Wigner fit to the data with E-^cs = —1.75795 nieC^ 
and r = 4.12 keV for the U-Cf system at i? = 20 fm. 

The density-of-states approach requires a fit of the background which is sensitive 
(especially with regards to the width) to the fit boundaries. 



3.3 Conclusion 

An important aspect of these methods is that they show how the mapped Fourier 
grid (MFG) method handles the supercritical resonance state which is important 
for the time- dependent work in chapter 6. The collision work requires the di- 
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Figure 3.2: The density of states Eq. (3.4) from a hermitian matrix diagonalization 
for the U-Cf {R = 20 fm) Hamiltonian obtained with a basis with = 3900 states 
and s = 6000, is shown by crosses (+). The hne shows a fit with Eq. (3.5), and 
yields E,^^ = -1.75794 meC^ and T = 4.17 keV. 

agonahzation of both subcritical and supercritical matrix representations of the 
hamiltonian. The above demonstrations show that the supercritical resonance is 
represented both in the density of states and in the wavefunctions of the eigen- 
states. Thus, no augmentation or alteration of the MFG method is needed in order 
to include the supercritical resonance state. 

While the calculations here used a large number of states, which would not 
be usable by the repeated diagonalization needed for the time-dependent work, it 
nonetheless shows that the eigenbasis does contain the supercritical resonance state. 
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4 Analytic continuation methods applied to 
supercritical states 



Analytic continuation methods take a real hermitian Hamiltonian and transform 
it into a non-hermitian complex Hamiltonian. A parameter governs this trans- 
formation, the analytic continuation parameter (ACP), (, and the original, real 
Hamiltonian is recovered when ( = 0. The goal of these methods is to attenuate 
the exponentially divergent resonance state, (p^csix), bringing it into the Hilbert 
space of the (analytically continued) Hamiltonian. The resonance state is not in 
the hermitian domain of the real Hamiltonian since it does not satisfy the boundary 
condition: (pres{x) — >• as x — > oo [26]. 

In order to obtain the best estimate for resonance position and width, the 
parameter (, should be chosen carefully since the resonance parameters will be 
C-dependent for a finite matrix representation. This is most often accomplished 
by changing ( and looking for stability in the resonance position £"1.08 and width 
F. These parameters, {Ej-es, T} can be plotted parametrically as ^-trajectories: 
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(-Ei.cs(C)) r(C)). The resonance parameters are determined from a non-hermitian 
eigenvalue problem with eigenenergy E'r = E'res ± zr/2 (positive sign for negative- 
energy resonance states and negative sign for positive-energy resonance states). 

Building a matrix representation of the analytically continued Hamiltonian, 
using the mapped Fourier grid method, at given values of ( and diagonalizing the 
resultant matrix yields a single spectrum of eigenvalues and a single eigenvalue 
for the resonance. By repeating this procedure for consecutive ^-values, a set of 
complex eigenvalues for the resonance is obtained. This set of eigenvalues is then 
used to find stable values of the resonance parameters as a function of (, yielding 
the best approximation to the resonance parameters (without extrapolation). Some 
of this work was originally published in [31]. 

4.1 Complex scaling 

Complex scaling introduces an analytic continuation of the Hamiltonian by a trans- 
formation of the reaction coordinate, 

r re'''. (4.1) 

It has been used extensively in atomic and molecular physics, and has been put on 
firm mathematical grounds by Reinhardt [27] and Moiseyev [26] and extended to the 
Dirac equation by Seba [32]. Recently it was applied to the case of Stark resonances 
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in the Dirac equation [33]. Increasing 9 in small steps beginning from zero results 
in a series of energy spectra where the continuum-state energies are increasingly 
rotated into the complex energy plane with the branch cut(s) as pivot points shown 
in Fig. 4.1 for the positive continuum (in the Dirac case the branch cuts are at 
imgC^). The bound states remain on the real axis and are not affected by the change 
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Figure 4.1: Illustration of the spectrum of a complex scaled non-relativistic hamil- 
tonian where the eigenvalues of the bound states ( x ) remain on the real axis but the 
eigenvalues of the continuum states (represented as a continuous line) rotate into the 
imaginary plane. The pivot point for the rotation is at the branch cut (continuum 
boundary). The continuum states rotate by 29 since the non-relativistic momentum 
is transformed by p — pe~*^ transforming E = p^/ (2me) — > e~^*^p^/ (2me). 



in 9. As 9 is increased from zero, the continuum states rotate into the imaginary 

plane by approximately 29 for non-relativistic states. This can be understood (non- 

relativistically) by the scaling in the momentum p = —i-^ — > pe^*^. 

The energy of the (non-relativistic) continuum states is approximately E = p'^/2 

(in natural units). When the momentum is scaled the energy becomes E — >• 
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g-2i0p2y'2_ jj^ ^]-^g same fashion, the relativistic energies rotate hj E ^ ±a/1 + p^e~^*' 
for small p [33]. The resonance states (located in the continuum) initially rotate 
with the neighboring continuum states as 6 is increased from zero, but then stabilize 
at the resonance energy (i.e., a.t E = E^^s ~ "^1^/2), and remain relatively stable for 
a range of 9 values (this is why the method is also referred to as the stabilization 
method). The closest approximation to the resonance energy occurs at the most 
stable energy eigenvalue along the ^-trajectory, which is found where | is mini- 
mized. Complex scaling (rotation) thus is a method to isolate the resonance states 
from the continuum. 

4.2 Smooth exterior scaling 

Smooth exterior scaling (SES) is an extension of the complex scaling (CS) method 
whose mathematical justification was developed, e.g., by W.P. Reinhardt [27], and 
N. Moiseyev [26]. In this thesis SES is extended to the relativistic Dirac equation 
for supercritical resonances. SES relies on the same justification as CS, but uses a 
general path in the complex plane that is continuous; when using non-continuous 
paths one refers to the method as exterior complex scaling (ECS) [26]. A simple 
path is obtained by rotating the reaction coordinate into the complex plane about 
some finite position, rg, instead of the origin. The transformation then has the 
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form, 



r 




(r — Ts) e*^ + Ts for Vg < r. 



r 



for r < 



(4.2) 



It offers the advantage of turning on the scahng at a distance which can be chosen 
appropriately for a given potential shape. It is natural to choose rg such that the 
"bound" part of the resonance state is not affected directly by the complex scaling 
(cf. section 2.3.2). 

This additional freedom introduces some comphcations not found in CS: unlike 
CS, SES does not always have a minimum in the |^^| curve as a function of 9, 
which makes it difficult to determine a stabilized 6 value, ^opt, for Eii{0). An ap- 
proximation can always be made by finding the cusp of the trajectory in {E^^s^ T} 
space [34], although this does not allow for a very precise determination of 6'opt 
compared with a minimum in |^^|- Alternatively, the parameters can be deter- 
mined from either the or ^ curves. This is because in practice it is usually 
found that at least one of them will have a minimum, with the best results for each 
parameter coming from its own derivative minimum, e.g. the minimum of ^ will 
determine the optimal value of 6 for T{6) but not for E^^^iO) [6]. Optimal results 
are obtained when the three derivatives have the same value of 6'opt, yielding the 
same values for i^ros and V. 
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4.3 Complex absorbing potentials 

The method of adding a complex absorbing potential (CAP) to the Hamiltonian 
has been used extensively in atomic and molecular physics [35, 36, 9, 10], and was 
put on firm mathematical grounds by Riss and Meyer [37]. In this thesis, it is 
extended to the relativistic Dirac equation for supercritical resonances. The true 
Hamiltonian is augmented by an imaginary potential, which makes the resulting 
Hamiltonian operator non-hermitian. In the case of the Dirac Hamiltonian a CAP 
is added ClS Sb scalar, 

HcAP = H-tV^W{r), (4.3) 

where ?7 is a small non-negative parameter determining the strength of the CAP 
and /3 is the standard Dirac matrix [21]. The function W{r) determines the shape 
of the absorbing potential and can be tailored for the problem at hand. Currently, 
the most common use of a CAP is as a stabilization method: one solves the system 
on an equally spaced mesh of values of rj and takes the minimum of 
the closest approximation to the true resonance parameters [37]. 

In the literature, one finds investigations of polynomial CAPs of the form 

W{r) = e(r - re) (r - rJ'' (4.4) 

where 6 is the Heaviside function and n is an integer. The CAP parameter 

allows for the turn-on of the absorbing potential outside of the "bound" part of the 
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as 

V=Vopt 



wavefunction [38]. 

Although CAP and SES appear as separate methods, they are not independent. 
It has been shown that a CAP can be transformed into CS [39] and that SES is 
related to CAP [8]. 

4.4 Results using analytic continuation methods 

Calculations were carried out using the potential in Eq. 3.2 for a U^^'^-Cf^^'^ quasi- 
molecular system for an internuclear separation oi R = 20fm where the atomic mass 
numbers Ay = 238 and Acf = 251. The results from CS, CAP and SES methods 
are compared here. For the present problem there is no barrier in the potential. 
However, as explained in section 2.3.3, one can find the region where the super- 
critical resonance state displays tunneling behavior by looking at the intersection 
of ii^res with the curve V{r)—meC^. For the present case, this occurs at ~ 2 (in 
units of h/m(,c). 

Due to the aforementioned difficulty (cf. section 4.2) in obtaining results using 
stabilization, smooth exterior scaling will not be compared to CAP or CS in the 
following sections. Its usefulness will be addressed in chapter 5. 
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4.4.1 Linear and quadratic CAPs 

When using a linear absorber (n = 1 in Eq. 4.4) it was found that the resonance 

parameters (-Eros; T) do depend on the choice of r^, but are least sensitive when 

Tc — ry- In fact, a minimum is observed for the width parameter, F, at = ry as 

seen in Fig. 4.2. For a quadratic CAP {n = 2), there is minimal variation in the 

extracted resonance parameters on the same scale. 
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Figure 4.2: The width F(1S(t) as a function of (in units of ^/mgc) for a linear 
(x) and a quadratic (+) CAP as obtained in a basis specified by A = 2000 and 
s = 6000 for the U-Cf (R = 20 fm) system. 



A few conclusions can be drawn from these results. The linear CAP, through 

the sensitivity of the result to the value of rc, is providing support to the notion 

of an effective potential barrier in the calculation at r ^ Tc- The quadratic CAP 
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calculation is deemed to be less sensitive to the value of due to the gradual turn- 
on of the absorber with continuous derivative. It leads to more accurate results and 
will be used in the subsequent comparisons. Both CAP calculations do provide, 
however, consistent results for the complex resonance energy position (not shown). 

4.4.2 Quadratic CAP and CS 

In Fig. 4.3, results for the different analytic continuation methods are compared 
for the U^^"'"-Cf^^^ system at i? = 20 fm. The resonance energy Ej-^g and the width 
r are shown as a function of the mapping parameter s (cf. Eq. A. 7), for a fixed 
mesh size of = 2500 collocation points for the CS method and A^ = 2000 points 
for the quadratic CAP method. The linear CAP results would be indistinguishable 
from the n = 2 results on the scale of the figures. 

It was found that the CS and CAP {n = 2) methods return consistent results 
within a parameter range of 500 < s < 8000 with small relative fluctuations. Both 
Ej-es and P are determined with a relative accuracy of better than 10~^. For values 
of the mapping parameter s outside of the ideal range the uncertainties increase, 
particularly so for the width. 

Also shown in Fig. 4.3 are two results for the same system from Ref. [1]. They 
deviate among themselves on the same scale, approximately, as from the present 
values derived from the CS and CAP methods. Not shown are the results from 
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Figure 4.3: The width r(lS(j) and mean energy Ej.cs{lSa) for the U-Cf {R = 20 
fm) system as a function of the mapping parameter s (cf. Eq. A. 7) with = 2500 
points for the CS method (x) and = 2000 points for the CAP method (+). The 
hues indicate values from Ref. [1], the phase-shift result is given as a solid line, and 
the truncated potential method result as a dashed line. 



the Breit-Wigner fits based on the diagonalization of the hermitean Hamiltonian 
matrix in the Fourier grid method (cf. section 3.2.3). We note, however, that E'res 
is predicted in close agreement with the CAP or CS data, while P is over-estimated. 
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For the projection method it happens to agree with the result from a phase-shift 
analysis performed in Ref. [1]. 

In order to demonstrate the convergence properties of the present CAP and CS 
results, Fig. 4.4 shows the dependence of F on the grid size A^. A fixed value of the 
mapping parameter [s = 1000) is used. The data indicate a systematic variation 
4.0849 



4.0848 
4.0847 
4.0846 
5:4.0845 

CD 

tn 4.0844 
4.0843 
4.0842 
4.0841 
4.0840 



■ + 




+ 


1 

+ 


1 

+ 


+ 

X 


X 

+ ■ 










X 












X 








_ X 




. X . 











1000 1200 1400 1600 1800 2000 2200 2400 2600 

N 



Figure 4.4: The width F(lScr) as a function of mesh size A^ using the CS (x), 
and quadratic CAP (+) methods with mapping parameter s = 1000 for the U-Cf 
{R = 20 fm) system. 



with A^ for the CS data at large A^. For A^ = 2500, it appears as if the width is 
established to four significant digits. The CS results fluctuate (on this fine scale) 
for A^ < 1750, before systematically approaching the CAP value. The CS result at 
A^ = 1250 agrees with the CAP result only by chance. The CAP results, on the 
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other hand, display better convergence at = 1000 aheady, with a deviation of 
less than 20 meV from the = 2000 result. The stability of -Eres) for different grid 
sizes, is on the order of 10~^ nieC^ for both methods meaning it is stable for all 
values shown. 

4.4.3 Other configurations 

These results are generalizable from the test case of a U-Cf quasi-molecule at i? = 20 
fm to other situations. The relationship between the resonance position and width 
has been explored. In Fig. 4.5, F is shown as a function of -Eres over a substantial 
range of system parameters. On the scale of the graph, differences between the 
present data and previous literature results cannot be noticed, i.e., they follow an 
almost universal curve. A number of the data points were obtained from different 
quasi-molecular systems at various internuclear separations with both point and 
extended nuclear models (in monopole approximation). Data from single-center 
calculations with one extended nucleus also follow the same relationship. 
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Figure 4.5: The width r(lScr) as a function of £'res(lScr) calculated by dif- 
ferent methods and using different supercritical potentials. The U-U and U-Cf 
data from Ref. [1] were obtained for the distances R = 16*, 16, 20, 25, 30 fm, while 
the single-center CAP/CS data from the present work were obtained for Z = 
195,193,190,187,184,181,178. The leftmost U-U and U-Cf points marked by R = 16* 
fm were calculated with point-like nuclei. 



4.5 Conclusion 

The quadratic CAP and CS are found to give results in close agreement for a 
wide range of parameters. All methods used gave very stable results as compared 
with the quoted results from Ref. [1] and the direct methods of chapter 3. The 
CAP method has been extended to the Dirac equation for the first time, and a 
quadratic CAP was found to be superior in determining the supercritical resonance 
parameters compared with CS. 
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5 Augmented analytic continuation methods 



Analytic continuation makes use of a parameter ( to turn the original Hamiltonian 
into a non-hermitian operator. This raises two problems: (i) one has to optimize 
( to find the best approximation to the complex resonance eigenenergy; (ii) one 
worries about the effect of the unphysical ( parameter on the final result. The 
optimization in step (i) is performed by stabilizing some measure (e.g., 
depends on the eigenenergy as a function of ( (cf. chapter 4). 

Recent work based upon the method of adding a complex absorbing potential 
(CAP) has demonstrated how to obtain more accurate results by extrapolating the 
complex £'/j((^)-trajectory to C = [12]. This idea is extended to the relativistic 
Dirac equation in this chapter using the complementary methods of smooth exterior 
scaling (SES) and CAP. The stability and parameter independence of the results 
is demonstrated. The extrapolation technique allows one to obtain highly accurate 
results for smaller basis size, A^; thus enabling the extension of the calculations 
beyond the monopole approximation to the two-center potential. Some of this 
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work was originally published in [40] and [41]. 



5.1 Pade approximant and extrapolation 

The goal of the analytic continuation methods is to make the resonance wave- 
function ipj-es a bounded function by choosing an analytic continuation parameter 
C > Ccrit- Although ipj-es IS an eigenfunction of the physical, hermitian, Hamiltonian, 
it is not in the Hilbert space, since it is exponentially divergent. For a sufficiently 
large critical value of ( = Ccrit (called 6'crit for CS, and r/crit for CAP), ipj-cs becomes 
a bounded function and is, therefore, in the Hilbert space of the physical Hamil- 
tonian [26]. Taking the lim^_»o -E'/j(C) always yields a real eigenvalue corresponding 
to H{( = 0)\E' = -E'_r(C = 0)\E' since H{( = 0) is hermitian when acting on bounded 



The authors of Ref. [12] proposed the following: instead of directly taking the 
C = limit, an extrapolation of a part of the trajectory, Eji{(), namely for ( > 
Ccrit, is used to obtain the complex- valued -Epadc(C = 0). The points used for the 
extrapolation are computed eigenvalues restricted to a region where ip^es remains 
in the Hilbert space (as represented by the finite basis). For the extrapolation the 
Fade approximant is given by, 



functions. 



-E'pade(C) 




(5.1) 
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where Pi and qj are complex coefficients [42] . It is useful to define Np = 2 [Ni + 1 ) as 
the number of points used in the approximant. The extrapolated value of -Epade(C = 
0) is given by po, and follows from a given set of ^-trajectory points, 

e, = {e, = Er{Q + (j - 1)AC), j = l...iVp} . (5.2) 

The smallest ^-value in the set of points used for the extrapolation is labeled Q and 
must satisfy Q > Ccrit- 

At ffist glance it would appear that new parameters have been introduced (A^p, 
(i) to replace (. However, it turns out that the extrapolation is very insensitive 
to Np for 4 < A^p < 12. While the extrapolation is also not very sensitive to the 
initial data point of the trajectory sample Q, we can optimize Q and the starting 
parameter (rc or r^) together and remove the effect of these two parameters on the 
final result. This is done by finding po as a function of Q, i.e. for different starting 
values of the extrapolation set e^. The distance between all the po's (at different C 
starting values) for different or re's are compared and a minimum in phase-space, 
{£'i.es,r}, is found. Averaging over the results for the two resonance parameters 
gives the best value and the error. This method allows for any number of or r^'s 
to be compared and even works when mixing SES and CAP results. 
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5.2 Smooth exterior scaling results using extrapolation 



In SES (as in the other methods) the energy Epi{6) forms a trajectory which de- 
pends on the calculation parameters {N,s,rs). Figure 5.1 shows the 6'-trajectory 
for = 250 (number of collocation points), s = 400 (mapping scaling parame- 
ter), A6 = 0.01 (choice of spacing for the trajectory) and = 2 (exterior scal- 
ing radius) and = 3 respectively as a sequence of squares and circles respec- 
tively. The values are displayed as the deviation from a reference value given by 
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Figure 5.1: The r(lS(T) resonance ^-trajectory from = 250, = 2 (□) and 
Ts = 3 (©) displayed as the difference from a reference calculation (cf. text) for 
the U-Cf system at i? = 20 fm in monopole approximation. The x's mark the 
resonance energy obtained by minimizing | | . The A^p = 8 data points used for 
the extrapolations are bracketed by the solid symbols, and a spacing of = 0.01 
was used. 
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£^res = -1.757930695(8)mec2 and T = 4.084798(12)keV, which was obtained by a 
large-size CAP calculation using = 3000, s = 400 and Tc = 2. The uncertainties 
were determined from calculations with different basis parameters, s. 

In Fig. 5.1, the Eii{9) results for small 9 values start in the upper right hand 
and abruptly change direction after a few points. This happens in the vicinity of 
the reference value given by the point (0,0). As 9 increases beyond the cusp in the 
trajectory, Eb,{9) changes more slowly, and the 6'-trajectory points become denser. 
Due to the finite number of collocation points A^, the actual value of ^crit (the 
minimum 9 value such that ipj-cs is bounded and properly represented within the 
finite basis) is higher than theoretically predicted for the N ^ oo representation 
[26], and is close to the cusp-like structure in the 6'-trajectory. 

Optimal, directly calculated, approximate values for the resonance energy are 
chosen by finding the minimum of ^ [26, 27], and are displayed as black 

crosses for both choices of r^. The solid lines represent Ap = 8 Pade extrapolation 
curves (cf. Eq. 5.1). The values of Si and Sn^ from Eq. 5.2, i.e., the bracketing 
points used to determine the Pade approximant, are indicated by filled symbols. 
The values of 6'i = 0.28 for the = 2 curve and 9i = 0.36 for the = 3 curve were 
chosen, because their -E'pade(C = 0) points were found to be closest to each other. 
Using other starting values for the Pade curves would result in indistinguishable 
extrapolation curves on the scale of Fig. 5.1 (as long as the starting point would 
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not be too close to the cusp in the 6'-trajectory). From the Fade curves one obtains 
the best approximation to the resonance parameters by selecting the end points 
(£^Padc(C = 0)), located close to the the origin (0,0), i.e., close to the = 3000 
reference calculation result. 

Although very different trajectories from those shown in Fig. 5.1 were obtained 
for different values, the Fade curves still extrapolated to nearly the same C = 
complex resonance energy (£'res+«F/2), as long as was chosen outside the "bound" 
part of the wavefunction and not too large, i.e. 1.5 < rs < 6. Larger values of rg 
pose difficulties for the finite-A^ computation as it cannot span the large-r region 
properly. 

Simple complex scaling (CS) is equivalent to SES with rg = and, thus, falls 
outside of the range of acceptable values. It yields results with a much larger 
deviation from the reference result, even when Fade extrapolation is applied. 

The scaling parameter, s, introduced by the mapping of the radial coordinate 
in the mapped Fourier grid method in Eq. A. 7, plays a dominant role as far as 
numerical accuracy is concerned. In Fig. 5.2 the magnitude of the relative er- 
ror in the width F(lScr) is shown as a function of s for SES calculations with 

= 500, while using two reference CAF calculations. The latter were obtained 
from the stabilization method with a large basis (A^ = 3000) and yielded val- 
ues of E,es = -1.757930695, F = 4.084798 for scaling parameter s = 400, and 
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E = 

-'-'TCS 



-1.757930702, T = 4.084808 for s = 750 respectively. Each SES data point 
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Figure 5.2: Magnitude of relative error for Fade extrapolated widths AF(1S(t) 
for = 500 as a function of the scaling parameter s, for the U-Cf system at 
R = 20fm in monopole approximation. The data were obtained using two reference 
calculations from the stabilized CAF method with = 3000 and = 2; red plus 
(+) symbols: s = 400; blue squares (□): s = 750 (cf. text). 

was obtained by finding the closest intersection of the Vg = 2, 3, 4 Fade curves from 

all starting 6's (for 6 > ^crit)? and taking the average of these three best values. 

The use of different values within the acceptable range would cause changes too 

small to see on this plot. The resonance position, E^^s, is much more stable with 

respect to s and therefore not shown. 

The technique of combining information from different Vg calculations in order 

to select the ideal subset of 6'-points helps with the following problem. For each 
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calculation with a fixed value of rg we looked at extrapolated values for the reso- 
nance width as a function of the start value 6i, and observed the deviation from 
the average value. It was found that numerical noise present in the data was min- 
imal for some range of Oi values (different for each calculation). This noise was 
investigated as a function of extrapolation order A'p. Combining information from 
calculations with different values of allows for the elimination of two sources of 
random uncertainties associated with finite-precision input to a sensitive extrapo- 
lation calculation. 

The results demonstrate that the extrapolated {9 = 0) SES = 500 calculations 
exceed the precision of the stabilized (non-extrapolated) = 3000 CAP calcula- 
tions. Table 5.1 gives the resonance parameter results (E^cs,^) for two different 
extrapolations, namely A'p = 4, 8 for basis size A^ = 500, and mapping parameter 
300 < s < 700, which corresponds to a subset of the data shown in Fig. 5.2. The 
standard deviation of the mean. 



where x is the average x-value, is given separately for resonance position and width. 
The results are encouraging since extrapolation provides a marked improvement 
over the stabilized CS method discussed in chapter 4. 




(5.3) 
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r (keV) 


a{T) 


8 


-1.75793073 


1.7x10^^ 


4.0848085 


1.4x10"^ 


4 


-1.75793072 


1.6x10-^ 


4.0848031 


1.4x10-6 



Table 5.1: Averaged resonance position and width from the extrapolated SES 
method with basis size = 500 using extrapolation orders Ap = 8 and Ap = 4. 
The values were averaged over the stable range 300 < s < 700. The standard 
deviation of the mean is given by a separately for position and width. 

5.2.1 Pade extrapolation of CAP trajectories 

Previously it was shown that a quadratic complex absorbing potential (ra = 2 in 
Eq. 4.4) gives better results than simple complex scaling (CS) when used as a 
stabilization method as done in chapter 4. Further investigation has found that 
the performance of SES represents an improvement over that of CS (provided a 
minimum can be found), and that it is competitive with CAP. 

Even for the same calculation parameters (A^, s, Tc), the ?7-trajectory for CAP is 
found to be different from the 6'-trajectory in SES. Figure 5.3 shows a trajectory 
pair for identical calculation parameters, namely A^ = 500, s = 600, and = 
Tc = 3 respectively. Even though the trajectories are rather different, the Pade 
curves extrapolate to rather close values (£'pade(?7 = 0) and £'pade(^ = 0)), i.e. the 
extrapolated resonance parameters for both methods are very close. When looking 
at plots of the A^ = 500 r] = extrapolation results for CAP as a function of s, as 
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Figure 5.3: SES 6'-trajectory with = 500, s = 600, rg = 3 (blue circles), and 
CAP //-trajectory for N = 500, s = 600, Vc = 3 (red diamonds) displayed as the 
difference from a reference calculation for the ISa state in the U-Cf system at 
i? = 20 fm in monopole approximation. The x's mark the approximate values 
obtained from the stabilization method. The reference value was obtained by Pade 
extrapolation of CAP calculations with A^ = 3000, A'^p = 4 for rc = 2,3,4 yielding 



-1.7579307062(5)mec2 and F = 4.08481303(8)keV. 



was done in Fig. 5.2 for SES, almost identical results are observed. In table 5.2, the 
average and standard deviation of the mean, a (cf. Eq. 5.3), for the CAP results for 
A^ = 500, 300 < s < 700 are given for two different Fade point numbers, A'p = 4, 8. 
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r (keV) 


(Tr 


8 


-1.75793072 


1.6x10^^ 


4.0848141 


3.1x10"^ 


4 


-1.75793072 


1.7x10^^ 


4.0848122 


5.4x10-6 



Table 5.2: Same as in Table 5.1, but for the extrapolated CAP method with basis 
size = 500. 

5.2.2 Comparison of CAP and SES results 

Both methods (CAP and SES) show similar behavior with respect to different 
aspects of the calculation. Both are relatively insensitive to the number of Pade 
points, A^p, used in the extrapolation. It was found that the results were unchanged, 
to the reported precision, for 4 < A^p < 12. The effects of changing the distance 
between the analytic continuation parameter (A^ for SES and A?7 for CAP) were 
similarly found to be small compared with s. 

It was found that one could choose optimal Ci start values for the Pade ap- 
proximation by minimizing the deviation between -E'r(C = 0) from different or 
Tc calculations. This is simple to implement, and works well for different values 
of A^p. In this way, one obtains results that are independent of the starting point 
for the rotation (rg in SES) or the imaginary potential (rc in CAP). The results 
are, therefore, very stable with respect to the analytic continuation parameters and 
depend only on the parameters from the mapped Fourier grid method. As shown 
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in Fig. 5.2 for SES (a very similar graph was obtained for CAP), there is a range of 
s for which results are stable. For larger basis size N the stable s-range increases 
making a judicious choice of s less important. 

We have averaged the results from both the SES and CAP methods over the 
stable s region, for a basis size of = 500 in tables 5.1 & 5.2. The standard devia- 
tion for the width (t{T) within either method is below 10~^, while the results differ 
at this level. It is, therefore, of interest to determine the reliability of the error es- 
timate which is based upon basis parameter variations using larger- A^ extrapolated 
calculations. Comparison of the width results with such an estimate based upon 
A^ = 3000 CAP and SES calculations indicates that SES and CAP converge to the 
same value, closest to the CAP value given in table 5.2. 

Concerning the most appropriate order for the extrapolations it is worth noting 
that the CAP Pade trajectories are rather straight in comparison with the ones for 
SES, and Ap = 4 might be more appropriate in this case. Nevertheless, no system- 
atic improvement was found when going to the lower-order approximation (which 
might be deemed more stable with respect to numerical noise in the trajectory 
points) . 
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5.3 Beyond the monopole approximation 

The ability to compute the resonance parameters to high precision with moderate 
basis size (e.g., = 500) allows for the inclusion of higher multipoles which are 
present in the two-center interaction. While the effect on the resonance position 
is expected to be small, the sensitivity of the width to computational details (cf. 
the different results discussed in chapter 4) indicates that some broadening of the 
resonance may occur when the two-center nature of the potential is included. 

To account for two-center potential effects the wavefunction is expanded using 
spinor spherical harmonics, Xk,^) 



(5.4) 



which are labeled by the relativistic angular quantum number k (analogous to / 
in non-relativistic quantum mechanics) and the magnetic quantum number /i [21]. 
The Dirac equation for the scaled radial functions, /(r) = rF{r) and g{r^ = rG{r), 
then becomes (^ = c = me = 1), 

^-^/. = -{E-l)g^+J2{x.,,\Vir,R)\x.,,)9. , (5.5) 

K=±l 

, ±CXD 

J^ + -g^ = (E + 1)/.- ^ (x-«,,in^,i?)|x-.,,)/. , (5.6) 
ar r ^—f 

ft=±i 

where V{r^ R) is the potential for two uniformly charged spheres displaced along 

the z-axis, which is expanded into Legendre polynomials according to V{r, R) = 
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Y^'^^Vii^r, R)Pi{cos 9) [21]. The monopole term, Vq, is given explicitly in Eq. 3.2, 
and for the present work we include the coupling terms required for the k = ±1, —2 
channels (i.e., the Vi{r, R){x±K,^,\Pl\x±K,^l) dipole and V2{r, R){x±K,^l\P2\x±ii,^,) 
quadrupole terms). In general, the k = 1, —2 channels (-P1/2 and P3/2 respectively) 
have the strongest coupling to the k = — 1 (S-states) and are therefore expected to 
have the largest impact on the supercritical ground state. 

Table 5.3 gives the IScr resonance parameters for different internuclear sep- 
arations for both the one- and three-channel CAP calculations in the uranium- 
californium system at an internuclear separation oi R = 20fm. The CAP calcula- 
tions were performed using a basis size of = 500 per channel with 300 < s < 700. 
The CAP parameters of Vc = 2,3,4 were used to obtain the most stable ?7-range 
for extrapolation, which was carried out with order A'p = 4. Values in parentheses 
indicate the standard deviation of the mean (cf. Eq. 5.3). It is clear that the accu- 
racy of the calculations - as indicated by the deviation - is much higher than what 
is required to measure the effect of the the P-state channels on the ISa resonance. 

The final column contains the relative difference of the width between the three- 
channel results and the monopole approximation. The correction due to the dipole 
and quadrupole potentials is seen to increase with internuclear separation R. This 
trend can be expected, as the overall interaction becomes less spherically symmetric 
in this limit. The effect is small, however, as the strongest contribution towards the 
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rLyllll ) 


Single-channel (k = — 1) 


Three-channel ( 


«: = -l,l,2) 


iri-ral 


Ercs (mc^) 


Ti (keV) 


Eres (mc^) 


Ta (keV) 


(xlO-2) 


16 


-2.00635363(3) 


8.148233(4) 


-2.00646180(3) 


8.150153(3) 
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18 


-1.87487669(4) 
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20 


-1.75793073(4) 


4.0848148(4) 
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4.087394(1) 


0.631 
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-1.65393272(3) 


2.708987(1) 


-1.65414448(3) 


2.711524(2) 


0.937 


24 


-1.5612122(1) 


1.6966694(3) 


-1.56144826(1) 


1.698947(1) 


1.34 


26 


-1.47820830(4) 


0.9877052(4) 


-1.47846358(4) 


0.989571(1) 


1.89 


28 


-1.40354329(5) 


0.5221621(4) 


-1.40381341(4) 


0.523544(1) 


2.65 


30 


-1.33603643(6) 


0.2420749(2) 


-1.33631779(6) 


0.242976(1) 


3.72 


32 


-1.27469286(5) 


0.0931429(1) 


-1.27498253(5) 


0.093639(5) 


5.32 


34 


-1.21868219(5) 


0.0271161(1) 


-1.21897787(4) 


0.027322(2) 


7.59 


36 


-1.16731491(6) 


0.0050374(1) 


-1.16761472(7) 


0.005091(2) 


10.7 


38 


-1.12001745(6) 


0.0004170(3) 


-1.12031983(6) 


0.000439(4) 


52.5 



Table 5.3: Averaged ISa resonance position and width from the extrapolated CAP 
method with basis size = 500, Ap = 4, for single-channel (monopole) and three- 
channel calculations as a function of separation, R, in the U-Cf system. The values 
were averaged over the stable range 300 < s < 700 using = 2,3,4. The digit in 
parentheses represents the standard deviation from the average and the last column 
gives the relative difference of the width obtained from three-channel versus single- 
channel calculations. 
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binding energy (and thus the supercriticahty) is provided by the monopole part. 
The dipole interaction contributions to the resonance parameters are a result of the 
relatively small charge asymmetry in the U^^^-Cf^^"'" potential. 

By further including couplings to the D-states (k = 2, —3 or D^i2 and -D5/2 
respectively), which includes an S — D quadrupole coupling, significant corrections 
to the supercritical resonance parameters are obtained. Table 5.4 gives the results 
of the five-channel calculation compared to the single- channel data. While coupling 
up to P5(cos(^)) are included, it is the quadrupole interaction between the 5" and 
D states which causes the change in the resonance parameters. Since the dipole 
interaction is small for the present system, the next-order effect after the monopole 
interaction is the quadupole coupling. The large changes to the resonance param- 
eters are due to the inclusion of five-channels resulting in the S'-channel coupling 
to the D-channel by the quadupole coupling (in three-channel the quadupole only 
coupled K = —2 with k = 1). 
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rLyllll ) 


Single-channel (/€ = —!) 


Five-channel (k 


= ±l,±2,-3) 


iri-rsi 


Ercs (mc^) 


Ti (keV) 


Eres (mc^) 


Ts (keV) 


(xlO-3) 


16 


-2.00635363(3) 


8.148233(4) 


-2.03582804(3) 


8.681022(8) 


65.38 


20 


-1.75793073(4) 


4.0848148(4) 


-1.79673776(4) 


4.650882(4) 


138.6 


24 


-1.5612122(1) 


1.6966694(3) 


-1.60546338(2) 


2.150177(9) 


267.3 


28 


-1.40354329(5) 


0.5221621(4) 


-1.45044923(6) 


0.794256(9) 


521.1 


32 


-1.27469286(5) 


0.0931429(1) 


-1.32257286(8) 


0.201075(4) 


1158.8 



Table 5.4: Averaged IScr resonance position and width from the extrapolated CAP 
method with basis size = 500, A^p = 4, for single- channel (monopole) and five- 
channel calculations as a function of separation, R, in the U-Cf system. The values 
were averaged over the stable range 300 < s < 700 using Tc = 2, 3, 4. The value in 
parenthesis represents the standard deviation from the average and the last column 
gives the relative difference of the width for the five- and single-channel calculations. 

5.4 Conclusion 

We can summarize the findings from this chapter as follows: the generalization of 
the extrapolation method from CAP to SES, and the combination of data from 
both methods allowed us to obtain highly accurate data for relatively small basis 
size A^ (utilizing results for an optimum range of the mapped Fourier grid scaling 
parameter s, as well as analytic continuation method parameters Tc, r^). With this 
breakthrough in computational efficiency we were then able to explore the sen- 
sitivity of the resonance parameters to two-center interaction effects. For larger 
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internuclear separations (near the onset of supercriticality) these effects are appre- 
ciable. The augmented methods are also applicable to all resonance phenomena to 
which analytic continuation methods are applied. 
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6 Time-dependent supercritical states 



While static supercritical resonance states due to a single nucleus are unlikely to 
exist (or be sufficiently long-lived if such a nucleus was created by fusion), it is pos- 
sible to gain access to this phenomenon using heayy-ion collisions near the Coulomb 
barrier. In these collisions, the atoms are taken to be fully-ionized for reasons men- 
tioned in chapter 2 and for the practical reason of theoretical simplicity due to the 
lack of electron-electron interactions (Coulomb or spin-orbit). 

A classical treatment of the nuclear motion by a trajectory R(t) generates an 
explicitly time-dependent Dirac Hamiltonian for the electron motion which does not 
conserve energy. Energy can be taken from the nuclear motion to excite electrons. 
In the Dirac many-particle interpretation it also means electron-positron pairs can 
be created (cf. section 2.4.1). The potential was previously given in Eq. 3.2 for 
each nucleus. The potential is modified by the substitution R — > R(t), since the 
nuclei are no longer static. 

Direct numerical integration of the Dirac equation over time as a partial differ- 
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ential equation in 3+1 dimensions is practically impossible with present computers, 
therefore another approach is followed. The method proposed in this thesis relies 
on the fact that the mapped Fourier grid method generates a complete set of eigen- 
states (spanning the subspace) for the time-independent Hamiltonian in an efficient 
manner. All equations are in natural units = mg = c = 1). 

6.1 Time-dependent Dirac equation with particle creation 

The supercritical resonance, which is the principle focus of this work, primarily 
influences the pair creation spectra for electrons and positrons. These spectra, and 
the total amount of pair creation, are, thus, the desired quantities. In contrast to 
Schrodinger quantum mechanics where one has unitary time evolution of individual 
orbitals, the Dirac theory is to be considered as a multi-particle theory. This is true 
even if the initial condition involves a single occupied (or unoccupied) state. This 
follows from treating the Dirac equation as part of a quantized field theory. It 
requires that the propagator is applied to all initial states which will then allow 
for the calculation of the quantities of interest, i.e. the observables. This enables 
the calculation of the electron and positron creation spectra which can then be 
compared with experimental data. 
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6.1.1 Quantities of interest 

The signature for the existence of a supercritical resonance state is to be found in 
the positron and electron energy spectra. This is due to the decay mode of the 
supercritical resonance state resulting in an electron in the ground state of one 
of the receding ions and a free positron with an energy close to the supercritical 
resonance energy (cf section 2.2.1). 

The probabilities of electron creation (rij with i > F) and positron creation (nj 
with i < F) in a state i are given by 



Here F is the Fermi energy (cf. section 2.1), \x\ ) and \x\ ) are positive- and 
negative-eigenstates of the single-center Hamiltonian respectively. The propagated 
state, |0^(t)), corresponds to \xu) at t — oo (as explained in appendix B following 
the derivation of Ref. [30]). Thus, Eq. 6.1 implies that the population of electron 
states (bound or free) is obtained from evolving initial states \xu) I'pyif)) fo^^ 
ly < F, i.e., for negative-energy states to times t > tf (when the collision is over), 
and projecting onto the electron states. Similarly, Eq. 6.2 implies that the positron 
population at energy E ^meC^ +p^/(2mc) is obtained by propagating the positive- 
energy initial states and projecting onto the positron continuum states. 




(6.1) 



(6.2) 
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For the positron spectra, the ground state (|Xis )) provides the dominant con- 
tribution of all initial positive-energy states, since it has the strongest interaction 
with the negative energy continuum (cf. section 2.4.1). The propagated ground 
state (|0is(t)) = f/(t, — oo)|xis'')) "where U{t2,ti) is the propagator) thus yields a 
first approximation to the complete positron spectrum (which has the contributions 
from all positive-energy states). This lowest-order spectrum is given by 



The initial IS state (|Xis )) is propagated through the collision. The resultant state 
(|0is)) is projected onto the single-center positron states ^))- One can then 
add the contributions from other states in the same manner. 

The first excited state can be expected to play a significant role in dynamical 
pair production (cf. section 2.4.1). It also plays some role in the spontaneous pair 
production despite remaining subcritical. This is caused by the fact that the time- 
varying potential couples the first excited state strongly to the ground state even 
for intermediate nuclear distances. 

Thus, a more complete positron spectrum is obtained by adding the the con- 
tribution of the initial first excited state (|X2S^)) the spectrum from Eq. 6.3. 
Similarly, the other positive-energy states may be added, one by one. This allows 
for an additional check on the calculation since the contribution from higher states 
should diminish. 




(6.3) 
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It is not possible to propagate all states above the Fermi level, and the problem 
is compounded by the difficulty to propagate free states, since, in principle, the 
spectrum is continuous. To complete the electron spectrum (Eq. 6.1) is particu- 
larly difficult, since all the states to be propagated are continuum states 
As explained in chapter 2, the low-lying continuum states {\E\ ^ meC^) are un- 
likely to play a role since their wavenumbers are relatively large, i.e., their coupling 
due to energetic collision dynamics is weak (small matrix elements). Furthermore, 
one should expect that states of very high energy, \E\ 3> meC^, should not play 
a significant role. These states represent scattering energies too high to couple 
into the changing potential on the basis of absorbing energy from the nuclear mo- 
tion. Therefore, only some range of intermediate-energy scattering states should 
contribute substantially. 

The propagation of continuum states is possible within the mapped Fourier grid 
method. This is due to the continua being discretized, with adjustable coverage. 
The density of states is determined by both the scale parameter s, and the basis 
size A^. The discretization, by a finite number of states, means that one deals with 
quasi-continuum states. 

These states do not represent a single energy, but a range of energies, as they 
are localized in a finite volume of space. The energy resolution is limited by the 
fact that quasi-continuum (wavepacket) states have a finite energy width. The 
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advantage of a discretized E > nieC^ and E < — nicC^ continua is provided by the 
ability to propagate the quasi-continuum states as initial states in exactly the same 
way as for the bound states. One has the freedom to provide adequate coverage for 
the necessary energy ranges in Eqs. 6.1 and 6.2. 

A weakness of the mapped Fourier grid method is that the resultant basis is not 
perfectly orthogonal (due to numerical noise), which necessitates changes to the 
expressions for the expansion coefficients (given by {Xi\(pu{tf)) for an orthogonal 
basis). However, this does not result in any change to Eqs. 6.1 and 6.2, since they 
are derived for the expansion coefficients (cf. appendix B). 

6.1.2 Approximating the propagator 

The time evolution of the state 4>u(t) during the collision is carried out by the 
propagator. 



(where T is the time-ordering symbol), which acts on a state at time to and generates 
the state at t. The propagation evolves the state according the the time-dependent 
collision hamiltonian, H{t). The direct application of the propagator is very diffi- 
cult. It is, therefore, approximated resulting in a two-step process whereby the state 
is propagated over a short time segment after projection onto a quasi- stationary 




(6.4) 
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eigenbasis. This is achieved by approximating the propagator integral, given by 

t 

center (6.5) 

for small time intervals as a sum of small, rectangular areas. 

The approximation proceeds as follows. Given an initial state of the system, 
|\l/(r,to)), the state at a later time is given by, 

|^(r, t)) = U{t, to)|*(r, to)) = T jexp (^-i ^ dt'H{t')^ | |*(r, to))- (6.6) 

To approximate Eq. 6.6, small time steps may be used so that H{t) is assumed 
constant and may be taken out of the integrand. The propagator for each time 
step then becomes 

[/(ti,to) =exp(-ziJ(ti)(ti-to)). (6.7) 

Constructing a basis for H{ti) allows for simple time propagation. This is 
carried out by the projection of the initial state, |\l/(r, to)), onto the basis defined 
by if(ti)|0m) = Em{ti)\(l)rn) ■ Then the propagator is applied to the </)m-basis. The 
resultant state is given by 

*(r,ti)) = t/(ti,to)|*(r,to)) =^amexp(-z^„,(ti-to))|0m). (6.8) 

m 

This can be visualized as approximating Eq. 6.5 by summing a number of small 

rectangular areas. In an eigenbasis, H{t) Ei{t) for the ith eigenstate. Figure 6.1 

illustrates an E{t) curve for the zth initial bound state. 
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Figure 6.1: Illustration of the propagator approximation in an eigenbasis. The 
(blue) curve shows the eigenenergy of an initial bound state as a function of time. 
The total left rectangle (red plus green) shows the left-point approximation to the 
integral in Eq. 6.5 while the bottom (red) rectangle shows the right-point approxi- 
mation. The rectangle with a triangle on the top shows the trapezoid approximation 
actually used in the current work. Translation along the horizontal (time) axis is 
executed by propagating the state; while movement along the vertical {Ei{t)) axis 
is the result of projecting the state onto the H(t + At) basis. 

As the nuclei approach each other, the eigenenergy of the ith state decreases and 
reaches a minimum when the nuclei are at closest approach. The total rectangle 
starting at to (green plus red) shows the left-point approximation used to approxi- 
mate the integral of Eq. 6.5. The right-point approximation (used in Eqs. 6.7 and 
6.8) is shown as the lower (red) rectangle. Ei{t) is symmetric resulting in both 
approximations being equal. 

The propagation of the state is visualized in Fig. 6.1 as a horizontal displace- 
ment. The projection of the state onto the H{t + At) basis is visualized as a vertical 
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displacement. Using the right-point approximation thus amounts to propagating 
first and then projecting, while the left-point approximation is the reverse. By 
summing all these rectangles an approximation to Eq. 6.5 is obtained. 

In the present work, a better approximation to the U (ti, to) operator is employed 
by taking the average eigenvalue of the previous and current time for Eq. 6.5 giving 



This approximation to a definite integral is known as the trapezoid rule and has an 
error of O(At^). It is illustrated in Fig 6.1 as the triangle plus the rectangle. This 
approximation allows for the use of a coarser time-mesh to save computational time 
and does not require much overhead. 

6.1.3 Difficulties associated with a non-orthogonal basis 

In appendix A.l the problem of an non-symmetric representation of the Hamil- 
tonian is discussed in terms of the time-independent problem. Time-dependent 
problems also require special attention due to the fact that the representation is 
non- symmetric, since this leads to a set of eigenvectors that can be non-orthogonal 
due to the numerical noise, but still span the iV-dimensional subspace (all eigen- 
values are unique and real). 

In order to accomplish the time propagation of section 6.1.2 the wavefunctions 
must be projected onto the new set of states, as shown in Eq. 6.8. The set of 




(6.9) 
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states spans the Hilbert space, i.e., no norm should be lost due to the change of 
representation. However, the expansion coefficients are not simply given by 
the inner product of the wavefunction and the ith basis state. For a state 

|^)=^a„|0„) (6.10) 

n 

the inner product coefficient is given by 

Crn = {(t)m\'^). (6.11) 

Inserting 6.10 into 6.11 the expansion coefficient, a„, is related to the inner product 
coefficient according to 

Cm = ^a„((0m|0„)) (6.12) 

n 

^ CLn = ^Cm{{(f)m\ 4>n))~^- (6.13) 

■m 

It is clear that the non-orthogonality must be taken into account using the inverse 
of the inner-product matrix. This adds some overhead to the computation, since 
the inner-product matrix must be computed and inverted, but this process is sig- 
nificantly faster than diagonalization. Therefore, the inversion does not affect the 
overall order of the computational problem {0{N^)). 

6.2 Bare uranium-uranium collisions 

There are many different collision systems which become supercritical at closest 

approach and could be used to demonstrate the existence of supercritical resonance 
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states. The shorter the decay time the larger the positron yield from the super- 
critical resonance, making it easier to detect. Therefore, the heaviest system which 
is accessible to experiments is the most desirable. Over the next decade the GSI- 
FAIR collaboration is planning to perform bare uranium-uranium merged-beam 
collisions. Thus, we concentrate on this collision system in the present thesis. The 
computational method can be easily adapted to other bare-nuclei collision systems. 

A calculation of the complete collision dynamics is far beyond the currently 
available computer resources, and thus other approximations are made to the sys- 
tem in addition to the propagator approximation discussed in section 6.1.2. The 
current work is limited to zero-impact-parameter trajectories. This allows one to 
employ cylindrical symmetry for the Dirac field, eliminating the need for magnetic 
subchannels in the matrix representation. For zero impact parameter the coupling 
matrix elements arise from a multipole expansion of the potential (cf. section 5.3). 

In the present work we restrict the potential to the monopole term. This de- 
couples all angular symmetries and we can solve the dominant S-channel {k = — 1) 
in isolation. These approximations break the relativistic covariance of the Dirac 
equation forcing one to make a judicious choice of reference frame. 
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6.2.1 Choosing a reference frame 

While aU states need to be propagated, the choice of reference frame has a dramatic 
effect on the results. The center-of-mass frame, which appears as an obvious choice, 
is not very well suited for the mapped Fourier grid method in the present work. 
This is caused by the representation of continuum states as a quasi-continuum with 
a non- uniform coordinate mesh (cf. appendix A). The low-lying continuum states, 
which become bound states as the nuclei approach each other, represent electrons 
with low wavenumber moving too slowly to react efficiently to the (relatively) fast 
collision. 

Therefore, in order for the calculation to be carried out in the center-of-mass 
frame the upper continuum must be well sampled and the states should be properly 
represented. This is due to the restriction to the monopole approximation of the 
two-center potential which does not localize states in the vicinity of the nuclei for 
large internuclear separations. The monopole approximation, in fact, is only useful 
at short distances, since the initial and final states are defined for the hamiltonian 
of a single nucleus. 

In contrast to the center-of-mass frame, the target frame of one of the nuclei 
allows for the treatment of the second nucleus as a perturbation. The successful so- 
lution of the hydrogenic problem (cf. appendix A and Ref. [13]) gives us confidence 
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that the initial states are represented properly. Thus, the initial states are taken as 
eigenstates of hydrogen-like uranium. By starting the calculation at a sufficiently 
large internuclear separation, -Rmitiaij the basis change from the initial hydrogen- 
like uranium states to the perturbed hydrogen-like uranium will not induce any 
transitions. 

There remains a subtle problem of the proper frame to use during supercriti- 
cality. In section 4.4.3 the relationship between the resonance energy and width 
is demonstrated. This is important since the eigenvalue is frame dependent: as a 
consequence the resonance width (i.e., the positron yield) depends on the choice 
of frame. The center-of-mass frame gives deeper resonance energies for the same 
configuration and is more appropriate as the second nucleus is too close to be con- 
sidered a perturber. It is also the frame in which the monopole approximation to 
the interaction is defined properly. Therefore, we make a change of basis at an in- 
ternuclear separation of Rc.h. = 0.1 (38.6fm) both before and after closest approach. 
This must take place before supercriticality sets in, which yields an optimal choice 
right before the system becomes supercritical in either frame {R ~ 36fm). 

6.2.2 The trajectory 

To demonstrate the supercritical resonance effect the trajectory should maximize 
the time the system is supercritical, as discussed in section 2.4.2. Therefore, a 
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Coulomb trajectory of moderate energy (740MeV total center-of-mass energy) is 
chosen. The mutual repulsion of the nuclei slows them down (Coulomb barrier) 
resulting in very slow nuclei near closest approach allowing for the system to expe- 
rience supercriticality. 

The trajectory is generated by solving Eq. 2.9 numerically using constant time 
steps within three distinct regions. When the nuclei are far away the time step, 
At, is increased by an integer factor. This decreases the number of steps needed 
in the far region and is justified by the slow change of the eigenstates at large 
distances. By examining the eigenvalue of the ISa as a function of the internuclear 
separation, it was found that for uranium-uranium collisions at 740MeV this region 
occurs around R > 8 (3089fm). The opposite is true when the nuclei are close. The 
time step. At, used is reduced for i? < 2 (772fm). The division of the trajectory into 
these three regions allows for an increased time resolution near closest approach, 
and saves time-intensive steps (diagonalization and projection) in the far region. 

6.2.3 Time scales 

At = 740MeV, the uranium nuclei will reach about 16.5fm (center to center) 
at closest approach. The supercritical resonance at closest approach will be the 
deepest and have the shortest life-time (cf. section 4.4.3). This resonance has a 
life-time of 392 zepto-seconds (10~^^s). The collision system will be supercritical 
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for 2.3 zepto-seconds, a time which is shorter by two orders of magnitude. 

Even though supercriticahty lasts only for a short time, the calculation must 
span a much longer time frame. The collision system takes 250 zepto-seconds to 
reach closest approach (when started from far enough away to not induce any 
transitions due to changing frames from single center to the target frame). After 
supercriticahty, the collision must be propagated for at least another 2 atto-seconds 
(10~^^s). This is because there are still transitions taking place. 

Given that the fastest-decaying supercritical resonance has a lifetime much 
longer than the collision time, one is interested in considering nuclear sticking, 
where the nuclei remain fixed at closest approach (cf. section 2.5). 

The sticking times used should be realistic, i.e., experimentally obtainable in 
deep inelastic collisions. Such times can be on the order of a zepto-second (a few 
zepto-second sticking has been detected [19] and theory predicts slightly more with 
fully ionized atoms [17, 28, 29]). A sticking time of a few zepto-seconds is still 
only a few percent of the resonance decay time which would, therefore, seem to 
escape detection. One would, therefore, not expect to see a sharp supercritical 
resonance peak superimposed on top of the broad dynamical positron spectrum 
(cf. section 2.4.1). 

The supercritical resonance creation spectrum, for sticking times (T) much less 
than its life-time, may be conjectured as a Breit-Wigner shape (cf. Eq. 1.1) with 
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r ^ h/T . The height of this peak should also be related to the sticking time, 
T, since longer sticking would yield more supercritical resonance decay. Dynamical 
pair production on the other hand, which is unchanged during sticking, will interfere 
with the resonant pair production. This interference may provide a suitable method 
of detection. 

6.2.4 Results for a fully-ionized uranium-uranium collision 

The collision of fully-ionized uranium nuclei following a Coulomb trajectory at a 
center-of-mass energy of 740MeV was computed. The parameter space of the work 
is large, but each parameter has been minimized separately reducing numerical 
uncertainty to below 0.5%. The center-of-mass energy of 740MeV corresponds 
to a closest-approach distance of 16.5fm. This results in a deep resonance with 
Ey:cs = — 1.56meC^ and T = 1.68keV calculated using the technique of chapter 5. 
The nuclei are estimated to have a radius of = 7.44fm (using /2n = 1.2 X 
with A =238 for uranium). The nuclei are then 1.6fm away from touching. All 
results are in natural units [h = = c = 1) resulting in an energy of E = ±1 
corresponding to E = ±meC^. 

Figure 6.2 display the differential probability, dP/dE, for the positron spectrum 
using the parameters: N = 1024, s = 700, Rq = 25 (9.65pm), Rf = 200 (77.2pm) 
and At = 1.88 (2.42x 10~^^s). In the far region (R > 8) we used 4 x At and in the 
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near region {R < 2) we used At/100. The (green) curve with +'s is the spectrum 
obtained by propagating an initial IS state using Eq. 6.3. Each of the +'s indicate a 
positron quasi-continuum state {\Xi ^)) used for calculating (nj) from Eq. 6.3. The 
(blue) curve with the x 's is also the positron spectrum, but following propagation 
of all bound states (Eq. 6.2 with the sum terminating at the onset of the positive 
continuum). The (red) curve approximates the results from Miiller et al. [3] shown 
in Fig. 1.3. Propagating all bound states adds less than six percent to AP/AE at 




Positron energy [mc ] 

Figure 6.2: The positron spectrum for fully ionized U-U collisions at Ecm = 
740MeV. The +'s are from an initial IS state only (cf. Eq. 6.3) and the x's is from 
having propagated all bound states. The (red) curve approximates the results from 
Miiller et al. [3] . The (magenta) dotted curve is a (scaled) Breit-Wigner shape with 
the parameters of the supercritical resonance at closest approach [R = 16.5fm). 



104 



any energy. The addition of the excited bound states does not alter the spectrum 
significantly. As was expected, the positron states near the \E\ ^mgC^ boundary 
do not participate. 

The results obtained are in close agreement with Miiller et al.'s close-coupling 
calculation [3] . The only notable difference is on the lower energy side of the peak. 
The peak position is in close agreement for both calculations. 

The (magenta) dotted curve is a plot of the supercritical resonance structure 
(a Breit-Wigner using Eq. 1.1 but scaled to the peak of the positron spectrum) 
at closest approach {R = 16.5fm). The dynamically calculated positron peak is 
very broad compared to the supercritical resonance and is not centered at its peak. 
There is no indication of the supercritical resonance in this spectrum. One may 
adopt the view that the positron spectrum is too heavily dominated by dynamical 
positron production for the resonance signal to be observed. 

While the energy density near the supercritical resonance position is high, it 
would not necessarily be enough to show the sharp resonance peak. This may lead 
one to speculate that the peak is simply unresolved. The results from chapter 3 
address this. Despite it being unlikely that an eigenstate will fall at the resonance 
energy {XE=Ercs^^ surrounding eigenstates together represent the supercritical 
resonance state. This is due to the continuum states of the mapped Fourier grid 
method being quasi-continuum states with a finite energy width. Thus, if a sharp 
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resonance peak was to be found in the spectrum the neighboring eigenstates would 
have shown it. 

6.2.5 Enhancement due to nuclear sticking 

Calculations for uranium- uranium collisions with nuclear sticking have also been 
performed. The same parameters as in section 6.2.4 have been used. The only 
modification is that the nuclei are held fixed at closest approach for a time, T, 
given in zepto-seconds (10~^^s). 

As long as the nuclei stick together propagation occurs with a frozen H{ti) only. 
No projection is done since the basis is unchanged. This results in static occupation 
probabilities until the sticking ceases. The basis must, therefore, be large enough 
to allow for the created negative-energy wavepacket to propagate to larger radii 
without distortion. With a basis size of = 1024 this limits the computationally 
feasible sticking times to less than T ^ 150 zepto-seconds for uranium- uranium at 
= 740MeV. 

Figure 6.3 displays the positron spectrum, {ui), for trajectories with nuclear 
sticking times of T = 2, 5, 10 ordered from lowest to highest peak. All the bound 
states were used as initial states. The (red) curve, representing a sticking time 
of T = 2, has its peak shifted to the left as compared to the T = spectrum 
by approximately O.lSmgC^ (76.6keV). Longer sticking times result in higher and 
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Positron energy [mc^] 

Figure 6.3: The positron spectrum for fully ionized U-U collisions at = 740MeV 
for trajectories with nuclear sticking times of T = 2, 5, 10 zepto-seconds (10~^^s) 
ordered from lowest to highest main peak. Each curve was calculated using all the 
bound states. The (magenta) dotted curve is a (scaled) Breit-Wigner shape with 
the parameters of the supercritical resonance at closest approach [R = 16.5fm). 
The points with horizontal error bars of the same color represent the conjectured 
width {h/{T + To) with Tq = 1.15 x 10~^^s) of the resonant peak. 

narrower peaks. The (green) long-dashed curve, representing a sticking time of 
T = 5 has a secondary peak ai E = 3.5mcC^. The (blue) curve, representing the 
longest sticking time predicted by other theory groups [17, 28, 29], T = 10, shows an 
even higher and narrower peak with a width of 0.45meC^ (229keV). It was analyzed 
at Rf = 290 (instead of Rf = 200) to reduce the numerical noise. 

This curve also has two visible secondary peaks at E = 2.50meC^ and E = 
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3.23meC^. The second peak is higher than the primary peak of the T = calculation 
shown in figure 6.2. The results shown are in close agreement with those of Ref. [3] 
shown in Fig. 1.3. 

At half of the respective peak's height, a width estimate of the resonance decay 
spectrum fi/{T + Tq) is shown in the same color as the curve. Here To is half the 
time the system is supercritical without sticking, chosen to estimate the supercriti- 
cal decay without sticking. The comparison shows that the width of the dominant 
peak structures are not explained perfectly. There can be various reasons for this. 
The estimate ignores the interplay between dynamical and spontaneous positron 
production. Furthermore, the calculation includes only the contributions from the 
propagated initial bound states. The inclusion of initial discretized continuum 
states with energy En has been attempted, but with limited success. These calcu- 
lations do not converge reasonably with increasing i?„. This failure is attributed 
to problems with poor orthogonality of these states (both among themselves, and 
to a lesser degree with the bound states). Thus, the present calculation does not 
include contributions from direct free positron-electron pair production. 

The (magenta) dotted curve represents the (scaled) resonance. As the sticking 
time increases the supercritical resonance decay will play a more dominant role. 
Indeed, we see in Fig. 6.3 that the peak is becoming narrower and more centered 
near the supercritical resonance energy as the sticking increases. As the width of 
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the resonant pair production spectra (estimated at T ^ h/T) decreases with longer 
sticking times the interference peaks also shift to lower energy. This is due to the 
increased dominance of the spectrum by the resonant pair production. At longer 
sticking times (20 < T <^ ^/r) very little interference will remain as the dynamical 
pair production should be dwarfed by the resonant pair production. Currently we 
know of no way to separate out the resonant and dynamical contributions. 

One can interpret the results as follows: naively one would think that the spon- 
taneous decay results in positrons emitted at the resonance energy E'rcs within the 
natural width F. This would be expected for sticking times of the order of T ^ h/T. 
However, realistic sticking times are substantially shorter and result in a structure 
comparable to the dynamical positron spectrum. We speculate that the compli- 
cated spectra in Fig. 6.3 are a result of the two pair-creation processes (resonant 
and dynamical) interfering in a complicated way. 

6.3 Conclusion 

The collision of fully ionized uranium-uranium has been solved in the monopole 
approximation. The results obtained are consistent with previously published re- 
sults [3]. All bound states (obtained in the diagonalization of the single-center 
hamiltonian) have been used as initial states and the continuum was well covered. 

In collisions with no nuclear sticking, there is no evidence of the supercritical 

109 



resonance decay. The positron spectrum is shown to be a wide peak away from the 
resonance peak (at closest approach). 

The enhancement due to nuclear sticking can be significant due the the nar- 
rowing and shifting it causes in the main positron peak. This effect is large and 
also results in secondary peaks for a sticking time of T > 5. Thus, supercriti- 
cal resonance states may be detectable by their interference with dynamical pair 
production. 
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A The mapped Fourier grid method ^ 



The mapped Fourier grid method has been proposed for the non-relativistic Schrodinger 
problem and has been apphed successfully in the context of atomic and chemical 
physics [43, 44, 45]. We have extended the method to the Dirac equation, and have 
tested it in Ref. [13] for the exactly solvable problem of hydrogenic states. We 
show in this appendix how it yields a matrix representation for the Dirac-Coulomb 
problem. 

A.l The mapped Fourier grid method for the Dirac equa- 
tion 

The coupled radial Dirac equations for the upper and lower components G and F 
for a spherically symmetric potential, V{r), are commonly written as [21], [22] 

,dF (1 - K)hc ^ , /A 1N 

- ch- ^ — F = (E- V(r) - mc^)G, (A.l) 

dr r 

^^dG^il + K)hc^ = {E-V{r)+mc')F, (A.2) 
dr r 

^Some of the following was originally published in Ref. [13]. 
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where k is a quantum number analogous to the angular momentum quantum num- 
ber I in non-relativistic quantum mechanics. Using /(r) = rF{r), g[r) = rG{r) 
and the Coulomb potential for a charge Z, in natural units H = rrie = c = 1, the 



equations become 



^ - If 

dr r 



dg 



E 



Za 

r 



I] 9, 



K 



dr r 



E 



(A.3) 
(A.4) 



where a = /he is the fine-structure constant. It is useful to write Eqs. (A.3) and 



(A.4) in matrix notation 



t Za _ '\\ _ (A. _ !1\ 
\ T J \dr r J 



9{r) 



V 



fir) 



E 



9{r) 



I 



fir) 



(A.5) 



The mapping for the Schrodinger-Coulomb problem from Ref. [14] is first mod- 
ified by introducing an additional scaling parameter s giving, 



r{6) = s9 — Aarctan 



se_ 
1' 



(A.6) 



For the relativistic calculation we improve upon this map adding a second-order 
singularity at ^ = vr, giving 



r{e) 



s6 — A arctan ^ 



(A.7) 



This choice, with A = 4000, was motivated in Ref. [14] by looking at the coverage of 

the classical phase space for an optimal representation of about 20 bound hydrogenic 
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states. The introduction of the singularity at 6* = tt does not modify the behavior at 
small 9, and thus the choice oi A = 4000 is also good for the relativistic treatment 
of hydrogen. The added scale parameter s allows the tailoring of the coverage in 
phase space in order to achieve an optimal representation for more bound states 
or for a combination of bound and continuum parts of the spectrum. The latter 
is important for time- dependent applications where both excitation and ionization 
are present. 

The singularity at ^ = vr permits to map the semi- infinite r-range to ^ G [0, vr). 
This feature allows for incorporation of the correct boundary condition as r — oo, 
which leads to improved orthogonality of the computed eigenstates. 

In terms of the independent variable 6*, Eq. (A. 5) becomes, 



Za 1 

r{e) 



( 1 d I K 

\ \J(,e) de r(6») 



/ 1 d ^ 

I j(6») de r{e) 



— + 1 



v 



firm 



\ 
I 



E 



( \ 



V 



where J{6) is the Jacobian of the transformation from r to 6* given by 



s6 — A arctan ^ 
2 -, ^ 



(A.9) 



This Jacobian weakens the singularity in the second derivative at the origin (due 
to the antisymmetric extension into the — r region) in comparison with the original 
mapping given by Fattal et al. [14]. 

A finite sine series truncated at order N is chosen to represent each component 
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of the spinor, i.e., 



v 



g{r{d)) 

fim) 



j \ J2n=ibnsm{ne) ^ 



(A. 10) 



Eq. (A. 10) when inserted into Eq. (A. 8) yields an analytic expression for the deriva- 



tive, 



Y.l=i «m sin(m6') 



sm nfc 



Y.l=i ^(^m cos(m6 



\ 



(A.ll) 



In order to solve the problem numerically the 6 variable is discretized. The number 
of points used for the 6 mesh is set to A^, yielding a square (2A^ x 2A^) Hamil- 
tonian matrix. The spinor is now represented aX 9 = 9i by an eigenvector of the 
Hamiltonian matrix with the two components combined into a vector of length 2A^. 
The 6i points are evenly spaced by A6' = The furthest point on the mesh is 

given by r^ax = r{NA6), which is controlled by the scaling parameters. Note that 
Tmax increases with s in a non-linear fashion. Thus, for given A^, s determines the 
number of bound states and pseudo-Rydberg states that are obtained, as well as 
their respective accuracies. 

The expansion coefficients in Eq. (A. 10) are obtained from projection by the 
discrete sums: 



N 



ak 



N 



(A.12) 
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and 



N 



N- 



(A.13) 



The spinor (A. 10) is therefore written exphcitly as, 



Ivti Sf=i ^m=i sin(m^j) sin(m^,)^(% 
' Ef=iEliSinM,)sinM.)/(^,) 



y AT+l Z^j 

Introduction of the cosine-sine matrix 



(A.14) 



TV 



Dij = — mcos{m9i) sin(m6'j), 



m=l 



simphfies the discretized Eq. (A. 11) to 

^ im ii E7=i ELi Mme,) smim9)g{9,: 



'ELEliSinK)sinM)/(^,; 



y N+l \d9 Z^j 



( 



J 



i^E7=iA,/(%) y 



(A.15) 



(A.16) 



where the index i refers to the fact that the spinor derivative is evaluated at r 
Ti = r{9i). 

The discretized Hamihonian is now written in block-component form as: 



+ 



Dij kSi^j 



(A.17) 



Due to the derivative terms, the representation of the Hamiltonian in the Fourier- 
sine basis is non-symmetric. In practice, one should be concerned about the appear- 
ance of complex eigenvalues mostly due to numerical errors. Thus, it is desirable 
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to find alternate symmetric representations. In the non-relativistic 2 order differ- 
ential equation this can be achieved by re-mapping the eigenfunctions before the 
Fourier grid method is applied [44]. While this approach guarantees real eigenval- 
ues, it is also possible to lose efficiency as the introduced effective potential has a 
more complicated structure than the original one. 
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B Particle spectra from a non-orthogonal basis 



The mapped Fourier grid method yields a non-orthogonal complete set of states 
due to numerical noise. The derivation of the lepton/antilepton number operators 
is performed while taking into account the non-orthogonality, thereby correcting 
for it and yielding more accurate results. The derivation follows that of Wells et 
al. [30]. 

B.l Basis states 

In the target frame of a two-body collision, the limit oi t ^ ±00 results in a 
single-center potential for the target nucleus. It is, therefore, in this limit that 
we define the final particle states. The particles in the two-center basis (for finite 
times) are viewed as quasi-particles since they do not directly correspond to final- 
state electrons and positrons. In order to differentiate between particle and quasi- 
particle states, a subcritical QED ground state, |$o(t)), is chosen. It is taken to 
be the time-independent state of the bare, target nucleus, with the Hamiltonian in 
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natural units {h = c = = 1), 

HT = d-p + p + V{f). (B.l) 

For the particle basis the states \xi) are defined as solutions to the time- independent 
Dirac equation for a single nucleus 

(d ■ p + /3 + Kargctlr)) Ix^) = \x^) • (B.2) 

The basis corresponding to the quasi-particles \(f)i{t)) are solutions to the time- 
dependent Dirac equation for the full collision potential, the two-center modified 
Coulomb potential V2 center -R) , given by 

(a-p + /5 + \/2eenter(r,i?))|0.(t)) = . (B.3) 

Since the states obtained using the mapped Fourier grid method are complete but 
are not perfectly orthogonal (due to numerical noise [13]) the inner-product matrices 
for each are needed. For each basis they are defined as 

S^, = (X.IX,), (B.4) 
A, = {Mmit))- (B.5) 

Lastly, to distinguish between particle and quasi-particle states we use different 
symbols for their creation and annihilation operators. The operators are the 
creation operators for positrons and electrons respectively operating on the time- 
independent states of Eq. B.2. The operators are the creation operators 
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for quasiparticles and act on the states in Eq. B.3. The i index on the quasi- 
particle operator corresponds to the index for the time-dependent state \(pi{t)) when 
t — > — cxD which is proportional to \xi) making i a good quantum number in this 
hmit. 

B.2 Operator relations 

It is the time-evolved vacuum state, |<l>o(^))) which is calculated in practice allowing 
for the direct application of the quasi-particle operators. We are interested in the 
number of particles (not quasi-particles) making a relation between the two particle- 
operators types necessary. To do this we use two expansions for the field operator. 

The field equation for the field operator, in the two-center collision 

system is given by 

(a-p + /3 + V^2ccntcr(r,/2(t)))|^i(t)) =z^|^,(t)). (B.6) 

Firstly, the field operator is expanded into solutions to the time-dependent basis 
\(pj(t)) (with the usual switch for negative-energy states) of Eq. B.3, giving 

r<F s>F 

The symbol, F denotes the Fermi level which is placed above the negative-energy 

states (making positron states denoted by r < F) and below the ground state. This 

definition works only for subcritical systems, but this poses no problem since we are 
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only interested in the final particle content when the nuclei are separated. The state 
\(pr{t)) is the rth initial state in Eq. B.3, i.e. the IS, 2S, etc. Since limt^_oo <Pjit) 
Xj-e"*^-'*, {xi} satisfies the appropriate boundary conditions. Therefore we can also 
expand into the static target basis {Xi}: 

I^^W) = \Xr) e-''^^* + 5^a, |xs)e-^^^* . (B.8) 

r<F s>F 

Equating the two expansions in Eq. B.7 and B.8 we get 

^2^1 \Xr) e-^^'-* + J2a. \xs)e-'^^' = J^/^J |0p(t)) + . (B.9) 

r<F s>F p<F a>F 

Multiplying on the left by J^uixA Sj^l yields 



J2^^r{Xu\Sf^^l\Xr) 

r<F V 


^-iFrt 


+ 


^^as{Xu\S^^l 

s>F V 




(B.IO) 


P<F y 




+ 


a>F V 




(B.ll) 


'^Y^'Y^^^rSll{Xu\Xr) 

r<F V 




+ 


s>F V 


(B.12) 


P<F V 




+ 


a>F V 


^o{t)) 


(B.13) 


r<F V 




+ 


s>F V 


-iEst 


(B.14) 


P<F V 


Mt)) 


+ 


a>F V 


Mt)) 


(B.15) 


r<F 




+ 


s>F 




(B.16) 


P<F V 


Mt)) 


+ 


Y^'^f^.s^lkxuy 

a>F V 


Mt)) 


(B.17) 



We have now obtained the particle creation and annihilation operators in terms 
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of the quasi-particle operators giving 

bi = E,,.PlS^,l{Xu \Mt)) e^^^* + E^'^^M^X. \Mt)) e'"""' (B.18) 
h= Ep,./?p5,7i(x.|0pW)*e-^'=* +Y,^lS^,l{Xu\Mt)r e-^^'^' (B.19) 

cr,u 

4= E,,./5p5,-i(x.|0p(t))*e-^'=* +J2^^<^Sk,l{Xu\Mt)re~^^^' (B.20) 

(7,1' 

where a>F,p<F,k<F for 6|. and bk while /c > F for a^, and a^. 

We now simphfy the notation by defining the expansion coefficients in a non- 
orthogonal basis. In such a basis the expansion coefficients are not given by the 
inner products of the state and the basis states, (<y9n|^), but are derived as follows, 

1^) = 5^aj0„) (B.22) 

n 

Cm = (B.23) 

= 5^a„((0^| (B.24) 

n 

^fln = ((0m I 0n))~^ (B.25) 
m 

Therefore, we can rewrite the creation and annihilation relations using 

akAt)=Y.^k,l{Xn\Mt)) (B.26) 

n 

(note that a with a single index is an annihilation operator but with 2 indices it is 

an expansion coefficient) where u is the state which is propagated, for instance the 
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IS ground state. This simplifies the relation and puts it in terms of the expansion 
coefficients, ak,y{t) with the result 

4 = Ep<F /5^«fc,pe^''^* + E.>i. «<xafc,.e^^^* k < F (B.27) 
h= Ep<F + E.>F k<F (B.28) 

«fc= Ep<F/^Kpe''''* + Ea>Fa<x«fc,-e''''=* fc>F (B.30) 

It is noteworthy that the phase-factors are with respect to the energy of the fcth 
state and are not part of the sums. This implies that they do not contribute when 
evaluating number operators, e.g. ($o(^)|&I^fe|$o(^))- 

B.3 Observables 

We are interested in the number of particles that will be created after the collision. 
Therefore, the particle number expectation value needs to be evaluated. Using 
the relations in Eqs. B.27, B.28, B.29 and B.30 we can write the particle number 
operator in terms of the quasi-particle operators and operate on the time-evolved 
vacuum state |$o(^))- 

We work out the positron number as an example. The number of positrons 



122 



created in the kth state due to the initial states p and a are given by 



Uk) = (Mt)\blbk\Mt) 



JE^t 



\P<F 



J2 «<xafc,ae*^'=* j 

a>F J 



KP'<F 



-iEut 



(t'>F 



where it is understood that k < F. The only non-zero term is aa\ 



{rik) 



J2 «A^x,<x' (-^oW 

ct,<t'>F 



[B.31) 
[B.32) 



T>F 



For example, the number of positrons created in state (at t — oo) by propagating 
the IS only is 

(B.33) 



{nk) = \ak,is\'^. 



These are the same expressions obtained when an orthogonal basis is used. Indeed 
the electron number also works out to 



[B.34) 



CT<F 



where k > F. Therefore, we need only be concerned with how the expansion 



coefficients are formed in order to use the usual results for the particle observables. 



123 



Bibliography 



[1] Joachim Reinhardt, Berndt Miiller, and Walter Greiner. Theory of positron 
production in heavy-ion coUisions. Phys. Rev. A, 24(1):103-128, Jul 1981. 

[2] J. Schweppe, A. Gruppe, K. Bethge, H. Bokemeyer, T. Cowan, H. Folger, J. S. 
Greenberg, H. Grein, S. Ito, R. Schule, D. Schwalm, K. E. Stiebing, N. Traut- 
mann, P. Vincent, and M. Waldschmidt. Observation of a peak structure in 
positron spectra from u+cm collisions. Phys. Rev. Lett., 51(25):2261-2264, 
Dec 1983. 

[3] U. Miiller, T. de Reus, J. Reinhardt, B. Miiller, W. Greiner, and G. Soff. 
Positron production in crossed beams of bare uranium nuclei. Phys. Rev. A, 
37(5):1449-1455, Mar 1988. 

[4] P. Krekora, Q. Su, and R. Grobe. Transitions into the negative-energy dirac 
continuum. Phy. Rev. A, 70(5):054101, 2004. 



124 



[5] P. Krekora, Q. Su, and R. Grobe. Klein paradox in spatial and temporal 
resolution. Physical Review Letters, 92(4):040406, 2004. 

[6] P. R. Certain N. Moiseyev and F. Weinhold. Resonance properties of complex- 
rotated hamiltonians. Mol. Phys., 36(6):1613, 1978. 

[7] U V Riss and H-D Meyer. The transformative complex absorbing potential 
method: a bridge between complex absorbing potentials and smooth exte- 
rior scaling. Journal of Physics B: Atomic, Molecular and Optical Physics, 
31(10):2279-2304, 1998. 

[8] Hans O. Karlsson. Accurate resonances and effective absorption of flux using 
smooth exterior scaling. J. Chem. Phys., 109(21):9366, 1998. 

[9] Robin Santra and Lorenz S. Cederbaum. Complex absorbing potentials in the 
framework of electron propagator theory, i. general formalism. The Journal of 
Chemical Physics, 117(12):5511-5521, 2002. 

[10] Satyabrata Sahoo and Y K Ho. The complex absorbing potential method (cap) 
to study the stark effect in hydrogen and lithium. J. Phys. B, 33(28):2195, 
2000. 

[11] H. Masui and Y. K. Ho. Resonance states with the complex absorbing potential 
method. Phys. Rev. C, 65:054305, 2002. 

125 



[12] Roland Lefebvre, Milan Sindelka, and Nimrod Moiseyev. Resonance posi- 
tions and lifetimes for flexible complex absorbing potentials. Phy. Rev. A, 
72(5):052704, 2005. 

[13] E. Ackad and M. Horbatsch. Numerical solution of the dirac equation by a 
mapped fourier grid method. J. Phys. A, 38:3157, 2005. 

[14] Eyal Fattal, Roi Baer, and Ronnie Kosloff. Phase space approach for optimizing 
grid representations: The mapped fourier method. Phys. Rev. E, 53(1):1217- 
1227, Jan 1996. 

[15] T. Cowan, H. Backe, K. Bethge, H. Bokemeyer, H. Folger, J. S. Greenberg, 
K. Sakaguchi, D. Schwalm, J. Schweppe, K. E. Stiebing, and P. Vincent. Ob- 
servation of correlated narrow-peak structures in positron and electron spectra 
from superheavy collision systems. Phys. Rev. Lett., 56(5):444-447, Feb 1986. 

[16] T. Cowan, H. Backe, M. Begemann, K. Bethge, H. Bokemeyer, H. Folger, 
J. S. Greenberg, H. Grein, A. Gruppe, Y. Kido, M. Kliiver, D. Schwalm, 
J. Schweppe, K. E. Stiebing, N. Trautmann, and P. Vincent. Anoma- 
lous positron peaks from supercritical collision systems. Phys. Rev. Lett., 
54(16):1761-1764, Apr 1985. 

[17] Toshiki Maruyama, A. Bonasera, M. Papa, and S. Chiba. Formation and decay 

of super-heavy systems. Eur. Phys. J. A, 14:191, 2002. 

126 



[18] P. Salabura, H. Backe, K. Bethge, H. Bokemeyer, T.E. Cowan, H. Folger, J.S. 
Greenberg, K. Sakaguchi, D. Schwalm, J. Schweppe, and K.E. Stiebing. Corre- 
lated e"*" e~ peaks observed in heavy-ion collisions. Phys. Lett. B, 245(2):153- 
160, Aug 1990. 

[19] H. Backe, P. Senger, W. Bonin, E. Kankeleit, M. Kramer, R. Krieg, V. Metag, 
N. Trautmann, and J. B. Wilhelmy. Estimates of the nuclear time delay in 
dissipative u -|- u and u -|- cm collisions derived from the shape of positron and 
5-ray spectra. Phys. Rev. Lett, 50(23):1838-1841, Jun 1983. 

[20] R. Krieg, E. Boek, U. Gollerthan, E. Kankeleit, G. Klotz-Engmann, 
M. Kramer, U. Meyer, H. Oeschler, and P. Senger. Reaction dynamics stud- 
ied via positron and electron spectroscopy. Phys. Rev. C, 34(2):562-575, Aug 
1986. 

[21] Walter Greiner. Relativistic Quantum Mechanics: Wave Equations. Springer, 
3 edition, 2000. 

[22] J. J. Sakurai. Advanced Quantum Mechanics. Addison Wesley, 1984. 

[23] Carl D. Anderson. The positive electron. Phys. Rev., 43(6):491-494, Mar 1933. 

[24] Richard L. Liboff. Introductory Quantum Mechanics. Addison Wesley, 1984. 



127 



[25] J. Rafelski, L.P. Fulcher, and A. Klein. Fermions and bosons interacting with 
arbitrarily strong external fields. Phys. Rep., 38(5):227, 1978. 

[26] Nimrod Moiseyev. Quantum theory of resonances: calculating energies, widths 
and cross-sections by complex scaling. Phys. Rep., 302:211, 1998. 

[27] W. P. Reinhardt. Complex coordinates in the theory of atomic and molecular 
structure and dynamics. Ann. Rev. Phys. Chem., 33:223, 1982. 

[28] V. I. Zagrebaev, Yu. Ts. Oganessian, M. G. Itkis, and Walter Greiner. Su- 
perheavy nuclei and quasi-atoms produced in collisions of transuranium ions. 
Physical Review C (Nuclear Physics), 73(3):031602, 2006. 

[29] Valery Zagrebaev and Walter Greiner. Low-energy collisions of heavy nuclei: 
dynamics of sticking, mass transfer and fusion. Journal of Physics G: Nuclear 
and Particle Physics, 34(1): 1-25, 2007. 

[30] J. C. Wells, V. E. Oberacker, A. S. Umar, C. Bottcher, M. R. Strayer, J.-S. 
Wu, and G. Plunien. Nonperturbative electromagnetic lepton-pair production 
in peripheral relativistic heavy-ion collisions. Phys. Rev. A, 45(9):6296-6312, 
May 1992. 



128 



[31] Edward Ackad and Marko Horbatsch. Numerical calculation of supercritical 
dirac resonance parameters by analytic continuation methods. Phy. Rev. A, 
75(2):022508, 2007. 

[32] Petr Seba. The complex scaling method for dirac resonances. Letters in Math- 
ematical Physics, 16:51-59, 1988. 

[33] I. A. Ivanov and Y. K. Ho. Complex rotation method for the dirac hamiltonian. 
Phys. Rev. A, 69(2):023407, Feb 2004. 

[34] G D Doolen. A procedure for calculating resonance eigenvalues. Journal of 
Physics B: Atomic and Molecular Physics, 8(4):525-528, 1975. 

[35] M. Ingr, H-D. Meyer, and L. S. Cederbaum. Potential energy curve of the x^ 
(Ty + u resonance state of f^ computed by cap/ci. J. Phys. B, 32:L547, 1999. 

[36] Robin Santra, Robert W. Dunford, and Linda Young. Spin-orbit effect on 
strong-field ionization of krypton. Phys. Rev. A, 74(4):043403, 2006. 

[37] U. V. Riss and H-D. Meyer. Calculation of resonance energies and widths using 
the complex absorbing potential method. J. Phys. B, 26:4503, 1993. 

[38] Imke B. Miiller, Robin Santra, and Lorenz S. Cederbaum. Resonances and 

pseudoresonances in a potential with attractive coulomb tail: A study using 

analytic-continuation techniques. Int. J. Quant. Chem., 94:75, 2003. 

129 



[39] Robin Santra. Why complex absorbing potentials work: A discrete-variable- 
representation perspective. Phys. Rev. A, 74(3):034701, 2006. 

[40] Edward Ackad and Marko Horbatsch. Supercritical dirac resonance parameters 
from extrapolated analytic continuation methods. Phys. Rev. ^4, 76(2):022503, 
2007. 

[41] E Ackad and M Horbatsch. New calculations for heavy-ion collisions with 
super-critical fields. Journal of Physics: Conference Series, 88:012017 (8pp), 
2007. 

[42] L. Schlessinger. Use of analyticity in the calculation of nonrelativistic scattering 
amplitudes. Phys. Rev., 167(5):1411-1423, Mar 1968. 

[43] C. Clay Marston and Gabriel G. Balint-Kurti. The fourier grid hamiltonian 
method for bound state eigenvalues and eigenfunctions. The Journal of Chem- 
ical Physics, 91(6):3571-3576, 1989. 

[44] V. Kokoouline, O. Dulieu, R. Kosloff, and F. Masnou-Seeuws. Mapped fourier 
methods for long-range molecules: Application to perturbations in the rb2(0^) 
photoassociation spectrum. The Journal of Chemical Physics, 110(20):9865- 
9876, 1999. 



130 



K. Willner, O. Dulieu, and F. Masnou-Seeuws. Mapped grid methods for 
long-range molecules and cold collisions. The Journal of Chemical Physics, 
120(2):548-561, 2004. 



131 



-1 - 



■^m2 

00 




20t§b0''5 - 4000 



6000 



ski 



LL 



_I_L 



_L 



_LL 



8000 

_LU 



10000 



_L 



n- L 



-j^M'^u 14000 
_ij HI I -\- 



1 0e@ajlB 2&000§e- 



%rmm 



ifiiii)ini))si3i;.i;i).(i).(i).(i 



liSim 90063-06 1 00000 



AE^s (mc2 ) 



